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ABSTRACT 

The stress level in the milled skin at the root of 
a tapered, multicell swept wing is predicted by means of 
a matrix-force method. The solution is achieved by 
minimizing the internal strain energy, and the results 
compared with experimental tests. A single loading, 
consisting of a nose-up couple, is applied to each tip rib. 
The loading is transferred from the tip rib to an idealized 
structure by means of simple torsion theory. The idealized 
structure represents the inboard one-half of each semi-span. 

Results indicate that an accurate solution of the 
stress distribution in the actual wing can be achieved 
provided that the root boundary conditions are preserved. 
An extension of the analysis is Suggested in orcer to 
more closely define the maximum accuracy ienent in the 


particular matrix-force method of solution. 
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INTRODUCTION 

Highly swept wines of low thickness to chord ratios 
were first built with conventional stringer sheet methods, 
This introduced the problem of how to successfully predict 
and then to relieve the severe concentration of stress along 
the rear spar and aft root section skin. 

The stress concentration was reduced somewhat by shift- 
ing to multi-cell construction. Multi-cell construction as 
used herein describes wings having closely spaced spanwise 
spars or webs, and relatively few streamwise ribs. 

Wing analysis by customary beam bending and torsion 
methods was unable to accurately account for the root 
area stress levels unless modification factors were used. 
This was due primarily to the neglect of bending-torsion 
inter-action. However, as flight speeds increased and 
aspect ratios decreased, the inability of beam theory to 
properly define chordwise deformations caused even greater 
concern since thermo-aeroelastic problems replaced static 
stressing in order of importance. 

In order to properly account for bending-torsion 
inter-action, chordwise curvatures, shear deflections, the 
increasing use of large cut-outs, and very heavy milled 
skin, several alternative methods of analysis have been 
developed. Without considering the older relaxation 


techniques and their iterative solutions, each depends 





upon the use of a digital computer for solution. Each 
attempts to achieve a mathematical model which exhibits 
a properly deformed shape having absolute displacements 
comparable to the actual case. 

The most prevalent of these methods treat unomst wete 
ure as an assemblage of elastic components, which allows 
matrix formulation of the solution in terms of the various 
ee theorems. Argyris (Ref. 1) has showm that these 
energy theorems derive from the two funcamental principles 
of virtual displacements and virtual forces, both of which 
originated in the work of Maxwell, Mohr, and Engesser, 
Since both principles are independent of elastic laws or 
the structural material, their application to non-linear 
problems or to structures where initial thermal strains 
exist 1S guite easy. | 

The "matrix-force" method leads to an analysis in 
terms of forces as the unknowns, while the "matrix- 
displacement" method leads to an analysis in terms of 
displacements. Either method provides for exact satis- 
faction of equilibrium and compatibility in a structure 
mede™up of discrete elements. Simplification of structural 
modifications during the preliminary design stage is an 
additional advantage which accrues from the use of discrete 


elements. 
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Nach method is again categorized as to whether 
‘lumped-parameter” or "finite element" idealization of 
plane panels is used. "Lumved-parameter" idealization 
as used herein allows distinet properties of the structure, 
such as direct stress, shear, bending and torsion, to be 
separated or concentrated at discrete locations. The web 
or panel may be attached to the adjacent flanges at 
corners, mid-points of the panel, or continuously. "Finite- 
element" idealization prescribes plane panels of either 
triangular or rectangular shape, to be attached at corners 
or nodes of the structure. In either idealization, panel 
werping must be avoided, even by slightly revising the 
geometry of the structure, in order to achieve a rigorous 
solution. 

Although a complete solution, giving both forces and 
displacements, may be obtained from either the matrix-force 
or matrix-displacement method, it has been shown in Ref. 2 
that if only forces are of interest, the matrix-force 
method is the most efficient. Conversely, if only 
deflections are of interest, the matrix-deflection method 
is the most efficient. 

The degree of redundancy in the force method is less 
than for the displacement method. However, the stiffnesses 
of the elements, as used by the displacement approach, are 


Sasier to obtain than are the flexibilities of elements 








required by the force approach. The force method has 
the additional disadvantage of neslectins Poisson's 
effect for that vortion of the axial-load carrying skin 
which is lumced with the beam flanges. 

Argyris (Ref. 3) has shovm that while the force 
methed is nearly always more suitable for fuselage analysis, 
the displacement method may often be more convenient for 
the analysis of complex wines, since in the displacement 
method structural "node points" specify the matrix size. 
This generally results in smaller matrices than would be 
required to obtain similar accuracy by the force method. 

Current literature abounds with instances where the 
two methods, each formulated several ways, have been 
applied to laboratory specimens. This allows some study 
of parameter effects, but allows no true indication of the 
method avplicability to an actual wing which may possess 
extreme and sometimes abrupt variations in skin thickness, 
discontinuous spars, massive ribs, skin and wel cut-outs, 
and taper. 

It is the object of this report to compare predicted 
stress and strain levels near the root of an actual wing 
(which possesses all of the aforementioned analytical 
obstacles) with the stresses and strains obtained experi- 
mentally. The specimen to be analyzed was obtained from a 


Mach 2.5 all-weather fighter project which was canceled 








by the U. S. Navy before the structural test prosram was 
completed. The matrix-force method of analysis is used, 
and consists of a "lumpved-parameter" structural idealization 
as proposed by Wehle and Lansing (Ref. 4). The matrix 
formulation of the solution is similar to the scheme 
proposed by Lang and Bisplinehoff (Ref. 5) and organized 
in considerable detail in Ref, 6. This formulation allows 
efficient use of the digital computer and at the same time 
allows intermediate results to be printed out for study. 
This characteristic makes changes to the structural ideal- 
ization relatively simple to incorporate, and causes the 
effect of such a change to be readily apparent. Thus from 
the outset it was recognized that a follow-on analysis, 
using the same digital computer program, would allow 
investigation of the boundary conditions, and indicate 

the merit of the idealization with only limited additional 
effort. 

This work was conducted at the U. 5S. Naval Postgraduate 
School, Monterey, California, during the 1962-1963 school 
year. The authors are deeply indebted to Professor C. H. 
Kahr, who conceived the need for such an investigation. 
Professor Kahr, together with Associate Professor U. Haupt, 
rendered invaluable assistance and encouragement throughout 
the project. The authors express their gratitude to 


Aeronautical Laboratory Supervisor R. E&. McConnell and 








assistants R. A. Besel and R. O. Cunningham, whose 
cooperative efforts were so essential in the procurement, 
fabrication and installation of all experimental apparatus, 
Appreciation is also expressed to Chance Vought Aircraft, 


Inc,, for their cooperation in furnishines essential data. 





II THEORETICAL ANALYSIS 
li General 
The structure chosen for the analysis was the wing 

center section of an F8U-3, as shown in Fig. 1 below. This 
section consists of seven spars and one post beam per semi- 
span, upper and lower milled skins, two massive pivot ribs 
closely spaced about the aircraft center line, two inter- 
mediate ribs, and two heavy wing fold ribs located at the 


extremities of the center section. 


Pivot Ribs 






Post Beam 





Fig. 1 
Sketch of F&U-3 Wing Center Section 


The first requirement for the analysis was to study 
the pictures and drawings of the dis-assembled wing in 
oraer to establish a basis for simplifying assumptions. 
These simplifications would be necessary to enable 


representation of the structure by a mathematical model 





ae: re 


consisting only of axial-load carrying bars and plane, 
shear carrying, constant thickness panels. 

The first simplification was introduced by considering 
the pivot rib to be the root location. The pivot rib was 
thus assumed to be infinitely stiff and free of chordwise 
bending. This step was taken to simplify the mathematical 
model, reducing by 16 the number of components to be 
analyzed. This would undoubtedly reduce the accuracy of 
the solution. However, the effect should be predictable, 
and once the reduced solution is achieved, incorporation of 
the center section (pivot rib to aircraft center line) 
would require comparatively minor additional effort. 

A post beam, shown in Fig. 1, extended from the pivot 
rib ("root") to the intermediate rib in the actual wing. 
Since no shear web was involved, the post beam was removed 
from the idealized structure and its flance area distributed 
to the surrounding structure. The posts themselves were 
ipnored. The remaining structure was then represented by 
axially loaded bars and flat shear panels of constant 
matmekness, as shown in, Fig. 2. 

It. was then necessary to determine the cross-sectional 
areas of the bars and the thickness of the shear panels. 
This was done for the bars at the cross-sections located 
Seesbations y, = 25.095 anduy = 61.95, Pig, 3. A bar end 


erea included the existing flange area, and the total skin 
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‘GEOMETRY OF STREAMWISE CUTS 
AT yy = 25.095 AND y,, = 81.98 


Wing Station y, = 25.095 





' | Wing Station y,, = 81.98 » 


Cor.8 
12.576 






Xp 212.669 
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area perpendicular to the bar, taken from the center of 

the panel just ahead to the center of the panel just behind, 
as shown in Fig. 4. Spar web areas were not included. 

The concentration of axial load carryine skin at the flange 
locations has the effect of removing lateral contraction 


of the covers from much of the analysis ( [l= 0). 





Ww 81.98 + 


Distribution of Cover Skin to Peam Flansces 


Halon, Lb 


Difficulty was experienced in determining the "average" 
skin thickness of each of the Ih cover panels due to the 
span-wise variation of thickness. Chordwise variation in 
thickness was easier to account for since the aspect ratio 
of the individual panels was on the order of four. 
Representation of the multiply tapered skin by panels of 
average thickness will undoubtedly penalize the accuracy 


of the analysis. 
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The un-tapered bars of the intermediate rib were 
formed by adding to the existing chordwise flanges. The 
added area consisted of a skin area eoual to the lecal 
thickness times a plate width (equal to the panel width) 


distributed evenly on either side of the rib, as shown 


in F ijggaeer 





thickness 
used 






rib augmented by this skin 


Meyege 5 
Distribution of Cover Skin Area to Rib Bar Area 
The dimensions of the resulting elements of the 


idealized wing are listed in Table I. 


Ze Loading Method 

After establishing the idealized model, it was 
necessary to establish a means of load application. It 
must transfer a load identical to the tip-applied load on 
the actual wing to the idealized model at the intermediate 


rib. As shown in Fig. 6, a nose-up (the wing is mounted 
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Pic, 6 


Application of Couple to Wing Tip 
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inverted) couple perpendicular to the C.I.B. was applied 
to each wing tip. 

Considering only the instrumented wing, this couple 
must be converted to 12 cover shears, seven vertical web 
shears, and 14 axial flange loads on a streamwise cross- 
section (intermediate rib) at station y,, = 81.98, as shown 
in Fig. 7. This was accomplished by assuming simple torsion 
theory would be satisfactory inboard from the tip to 
approximately the intermediate rib. This assumption was 
based upon tests reported in Ref. 7. 

The shear flows induced at two cross-sections 
perpendicular to the C.I.B., and located at stations 


Yw = 74.3 and yw = 98.7, were determined from 


T= )2Aq (1) 


The calculations and results are shown in Appendix 
A. A linear interpolation, using the planform dimensions 
from Fig. 8, gave the shear flow for each cell, as oriented 
perpendicular to the C.I.B. These individual cell shear 
flows were then resolved in stresses along, and perpendicular 
to, the streamwise intermediate rib et station y,,-= 61.98. 
From Mohr's circle of stresses, shown in Fig. 9, there is: 


(90=5 20 e) 


nN fs 


f,= £, cos @ 


@© 
il 


Sin @ 


Fh 
tt 


15 





i= - => 's 
ees eS — — - - 
= w>- my 
' ' 
‘/ 
> 
=> 
| 
‘ ' ‘= 





FIG. 8 


CENTER SECTION OF F8U-3 PORT WING 
(Main Structural Elements) 
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Fig. 9 


Mohr's Circle of Stresses at Intermediate Rib 


Applying these stresses to the panel edge along the stream- 
wise rib, and taking average panel thickness t and edge 


distance ws; the forces F,, and F, due to streamvise shear 


Pp 
stress and stress normal to the rib, f, and fo» are: 


Fy 
F 


f, wt =f, sin © wt = ag w sin 6 


Lo wte=f, cos @ wt = q, w cos @ 


p 


It can be seen that the streamwise component may be 
carried by the rib flanges. However, the force component 
perpendicular to the rib, Fy, in Fig. 10, cannot be carried 


by the shear panel and must be separated into 4a streamwise 


ne/ 











Fig. 10 


Component of Shear Force 
Perpendicular to Intermediate Rib 


component, Aqn, which is additive to the force F,, and 
a component Ap taken axially by the adjacent flanges, as 
shown in Fig. ll. 


AP 


Nn 


Fig. 11 
Vectorial Relation Between Ad, and Ap 


One-half of theAp component was assigned to each adjacent 
flange. The various load components and the resulting 


applied loads on the idealized model are given in Table II. 


18 











Voy Pu 
JO MalA papojdx 7 


6! = 
d g/ “Oly 
ee a 
Sly ne 
Os ee ae eres 
= = iP 
Z a 
D 8219 / , 
DsO/t 7 
=, 
humane 























La 
e Sy 
Qei gs, Ry) 
® 
- 
8 
& = 
xy Qo 
ke) oO 
S) = 
“ 
NX 
~ te 
er bo 
$v 65/07 
——— 
w 
na J 
| n 
8 
~~ 
Fy t $n 
ee oe 
oe 
SS 
PU L/H6 
$£y — Fe 
£0 29/6 
; | 
vEY : fn 
SS ae 
'p0t28 


[9 











8p 95988 


mMaA papoj/ax7 


E/ 614 





2O 








ce Element Loads 

Loads were asSigned to the individual structural 
components as shown in Figs. 12 and 13. Only the redundcants 
and sufficient loads to define the individual member 
flexibilities and general state of stress were included. 

The statically determinate system was defined by 
making cuts inboard of the intermediate rib at station 
yw = 861.98. By "cutting" the five interior beam webs 
and all flanges except the lower front and two reer flanges, 
a minimum determinate and stable structure was obtained. 
A cross-section of the determinate structure is shown 


i rag, lh. 
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Fig. 14 
Cross-section of Determinate System 


21 





SS ——- => —_, 
i 4 eel” 





Since the intermediate "tip" rib at station y,, = 81.98 was 
considered to be determinate, there were 16 redundant and 
33 determinate loads. The load numbering system shown in 
Fie. 13 was established in order to facilitate organization 
of the analysis method, and to assist in conditioning a 
matrix for later inversion. Applied loads are therefore 
numbered P, through P33, and the determinate reactions are 
numbered aj through 233: The redundant loads were then 
assigned numbers a3, through ang: 

Shear panel edge loads were proportioned by first 
assigning a unit shear force to the outer edges of all 
cover panels and beam webs. Unit shear loads wene 
assigned to the forward edge of the rib webs. By 
representing a typical panel as shown in Fig. 15, where a; 


is the unit load, there is: 
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Fig. 15 


Proportioning of Panel Edge Loads 


The geometry of the panels and the ratios necessary 


for determining the loads are shown in Table III. 


be Bending Torsion Interaction 


Interaction between bending stresses and cover shears 
was accounted for by the method reported in Ref. 4. Swept 
panels were assumed to have parallel edges, as shown in 
Fig, 16. A rectangular panel was then formed with a constant 
shear flow equal to the average of the end shear flows, as 


given in Fig. 13. 
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Allowance For Interaction of Bending Stresses and Cover Shears 


The triangular segment remaining wes placed in 
equilibrium by reacting one-half its spanwise comvonent 


at each spar cap by the load P, where 
P = qw sin/\p 


Interior spar caps receive a contribution from panels on 
either side. These loads are shown by dashed arrows in 


ie, 13. 
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oe Equilibrium Equations 


Since the redundant loads were to be obtained by 
minimizing the internal strain energy, 33 independent 
equilibrium equations were required to establish the load 
distribution in the determinate structure. Because there 
were more than 33 members in the structure, some choice 
was available in writing the equilibrium equations. It 
was decided to avoid where possible the writing of equations 
for those members or combinations of members where it was 
likely a small load would be applied or resisted. 

The first six equations were obtained by expressing 
the equilibrium of the entire model in relation to the 


axis system chosen: 


ee eree= 0 
i) 22 
3.) LF = 0 
4.) LMxx = 0 
a a =o 
6.) Me, = 0 


Using the right-hand rule, each applied load and its "external" 


reaction may contribute incremental moments according to 


AMzz 


2p 








Cover sheers at either end were assumed to act at the 
forward or top edge of the panel or beam web. 

The delta moments are then summed and set eoual to 
BELO. 

The next eleven equations were obtained from axial 
equilibrium of the 11 redundant beam flanges, elements 16, 
18, 20, 22, 2h, 29, 31, 33, 35, 37, and 39. 

Ten more equations were obtained from the rib at 
station yw = 81.98 by considering the eouvuilibrium of inter- 
beam segments between the ReB. andwR. I. Ba, at. eB.-A, 1. Bie 
ile BweC.L.B., C.L1.B.-FW.1l.B.,eand MWD Beebe. Beeeor a 
rib segment, looking inboard from the tip, the moments about 
the y-axis (lower forward corner of the segment) and forces 
in the x-direction were set equal to zero, giving two 
equations per segment. 

Following Ref. 5, the vertical shear loads in the beam 
webs, applied at station y, = 81.98 by loads Pj through Py, 
were assumed to act through seven imaginary posts of unit 
cross-section, as shown in Fig. 12. The final six equations 
were obtained by setting all but the F.L2. post in equilibrium. 

This gave a total of 33 equilibrium equations. Each 
equation was arranged internally so that only the non- 
redundant (aj, through 433) loads appeared on the left side 
of the equation. The applied loads (P, through P33) and 


redundants (a3) - 4,9) appeared on the right side, 
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6; Matrix Formulation 
rollowing the matrix formulation given by References 
12 and 13, all at the equilibrium equations can be expressed 


by 


}] 33x33 tenth ; Pl as CE} Ox a 


Solving equation (1) for the non-redundant loads a,,. 
requires the inversion of the matrix of coefficients, [A] ° 
Therefore the row order of the 33 equilibrium equations was 
arransed to favorably condition the [a] matrix by either 
olacing the large terms on the main diagonal Or symmetrically 
placing groups of terms about the main diagonal. This in 


turn fixed the arrangement of the [B | matrix. The [A] and 


the [B| matrices are shown in Tables Iv and V. 


Since there are 16 redundants, 16 additional equations 
must be obtained by minimizing the strain energy of the 


systen. 


From Refs. 4 and 5, the internal strain energy of the 


structure may be written as 


U = = [a4 un ay, | F {24 ee an} (2) 


ay 








where clements of the symmetric IF matrix are the sums of 
the various member flexibilities. These flexibilities were 
obtained as shown in Appendix B. The [F| matrix was 
designated luc for the digital computer program, and 
BmoreLner with the B and [3] matrices, formed the ,basic 
structural input to the digital computer program shown in 
Appendix C, 

Referring again to equations (1) and (2) it is necessary 
to express the strain energy in terms of the external forces 
and redundant loads. 

Equation (1) was re-arranged by taking the inverse of 
matrix [A] , and post multiplying the inverse la~2| by 


matrix [B | » giving the equation 


{nx} 35 = [av lal Patel” A lasaug ‘| (3) 


r 


The column matrix {a art and the ca] matrix were 
increased to 49xl and 49x49 dimensions respectively, by 


simply adding the trivial equations 


a = a 
oe ol 
aug = alg 


The augmented matrix [ea] was then redesignated the 


g matrix, 
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The transpose of matrix [a] was then taken, allowing 


the strain enercvweraicx lH to be formed: 


Lien : (° _— F 49x49 G 49x49 ies 


By substituting equation (5) into equation (2), the internal 


strain energy may now be written as 


U = Py oe P32. a3h++830], [H] {P,. -P 234+ ayo} (6) 


The strain energy is then minimized by differentiating 
equation (6) and setting the result equal to zero, according 


to 
= il) (7) 


Since [x] is a symmetrical matrix, it may be shown that 
the operations of equation (7), as performed upon the expand- 


ed form of equation (6), may be expressed simply by 


— -1 
Oe.. P “te a 1? 33x11 ie 


The sub-matrices [Ho and Hoa | come from the 


partitioning of [H | according to: 


aS, 








applied load= 33 


1) 
al 

So 
ue) 
A 

=: 

4) 
i 


r = redundant=16 (9) 


By defining the producto. [Ho | x [Hoy | as: 


» a 
BR a3 ee | aul aad — 


tie answers to the redundant loads, as given in the computer 


program of Appendix C, are: 


Ph (8b) 


a Says : EH axa 


The total loading system is then determined from 


m2 ‘ se Pl a 
where the matrix [S| is constructed from [EN] and iG] by 
Be 33 
Is} = [¢] ~ eee (10) 
LOxh9 [EN] 16x33 


49x33 


Finally the internal strain energy is written in terms © 
of the applied loads and the flexibility influence 


coefficients ic as: 
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IP, eee (11) 


U=4 [Py ...-P3/[c | 


33x33 331 


By substituting equation (8b) into equation (6), the 


matrix of flexibility influence coefficients is given by 


fo] = 


Hy) - [Ryo] fen] (12) 


im 


ae Computer Programming 
The solution of equations (1) through (12) was adapted 


by a FORTRAN program to the CDC 1604 computer. The variable 
names and abbreviations, flow charts, programs and sub- 
routines are given in Appendix C. 

Considerable effort was required to establish the 
desired program within the basic computer storage capacity. 
Any future additions to the main program would require either 
use of machine language or an input-output program on 
peripheral equipment. 

Print-outs throughout the program were made in order 
to establish a means for continuing accuracy checks. 
Although a complete hand solution of one set of data is 
possible in many computer programs where merely repetitious 


iterations are done, here only equilibrium of elements and 
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standard machine inversion checks could be made. 

The first machine check was performed by writing a 
short auxiliary program namec CHECK, This confirmed the 
accuracy of the inversion of the matrix [A] as accomplished 
by the single-precision Gauss 3 subroutine. The matrix 


[cxmcx] was formed from 


wy TP ie 


[ CHECK] _ * 


= 4A 
B)oP 5/8 33x33 
33x33 
and is shown in Appendix C. 

The second check was made by comparing the member 
loads of the statically determinate structure With values 
irom the [ca | matrix shown in Table Vif. For example 
consider the statically determinate structure shown in 


eae. 17. 





aben AS 7/ 


Loaded Determinate Structure For ca Check 
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From statics a value of Ao, may be obtained for a unit 
load P97, This should compare with the value CJ 55 7 in 
3 


the expression 


P 
{Snz| 33x1 Gel 33x49 =| Lox ae 


where all terms in the column load matrix are zero except 
7. 
Setting the y-component of ajo equal to (aza)y, by 


statics: 


(azo)y = Py x (81.98 - 25.095) = 7.900694 Ph 
7. 20a 


209114 


From equation (14), and the matrix [cu] 


Ann = Py (8,8660330) + «+. + Py (8,8660330) 


‘al 


which checks well to the fourth decimal. 
By similar reasoning, the three non-redundant flanges, 


elements 14, 26 and 27, should react directly the loads 


Pon, Pobs Pa36 


Applied loads Reaction by ca] matrix 
Pow = 1.0 €o7 = 1.0000558 
Pog = 1.0 oo = 1.0000268 
P33 = 1.0 a3 = 1.0000238 
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Again considering a unit load at Poy the total shear 
load should be carried by the rear beam web, ag. From the 


ou matrix, a unit load at Pr gives 
ag = 1,0103419 


Consequently, the cover panel shear loads should be 


zero. From the [ea | matrix, 


aq = .0000646 ay5 = 0000437 
219 = 0000592 aig = .0000471 
11 = .0000550 ay7 = .0000511 
212 = .0000510 43 = .0000550 
a43 = .0000471 ayg = .0000591 
ay), = -0000438 Ang = 0000645 


The next check was to determine the state of 
equilibrium of each flange and post. For example, the 
front post of the intermediate rib, station y, = 61.98, 


mice, 12, gives 


1.C000000 = .5981273 - «4016727 = .0000000 


for a unit load at Pi: For the seventh post, the result 
from a unit load at Py is also zero. A maximum deviation 


from true equilibrium of 7.9 percent existed at the second 
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post from the front. It should be noted that no 
equilibrium equation was written for this post. 

Referring to Fig. 13, a unit redundant load was applied 
in turn to each redundant flange. The greatest deviation 


from equilibrium was in element 39, where 


Error = = 0.0C000043 1b 
[ca] 


Table VII contains the results of the CJ matrix 
equilibrium checks, which completes the check of the basic 
equations and computer operations on them, as well as the 
validity of the matrix inversion sub-routine. 

A print-out statement was used to determine that proper 
partitioning of the [HI matrix occured. Since computer 
storage space was not available, hand calculations were 
performed to confirm the proper inversion of Hp» | ° 

An interim calculation according to equation (8b) 
gave the redundant load matrix {AR} . One of these values 
was computed by hand as a random check and found to be 
correct. All 16 values were then checked with the final 
16 loads of the {AT} matrix obtained from equation (9). 

This completed the computer program checks. The 
program was then run twice and answers compared in order to 
detect any computer core errors. There were none. 

The flexibility influence coefficient, stress influence 
coefficient, and load matrices, ic|, [S| , and {AT} , are 
listed in Tables VIII, IX, and X. 
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§. Theoretical Results 

Compared with ordinary beam theory solutions, the 
accurate determination of stresses in an idealized structure 
by a matrix force method requires considerable engineering 
judgment. It has been pointed out by Warren (Ref. 9) that 
errors of one-hundred percent may occur in interpreting the 
results of such a solution. Therefore the original 
formulation assumptions were used in determining the results 
of this analysis. 

The investigation was narrowed to determining the 
average panel stress at station yy = 34.5 from the loads of 
the {as } matrix. For each panel, the maximum normal stress, 
the angle of the principal axis, and the maximum shearing 
stress were determined. From the components shown in Fig. 
18, and noting that the direction of shear flow reflects 


the minus sign indicated by load ag in the {AT} matrix, 





Fig. 18 


Theoretical Root Panel Stresses 
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the axial loéd components of the skin were determined. 
Using the solution formulation data of Table I, the bar 
loads (excepting a,,.) were re-apportioned by area ratios to 
the panels from which the skin came. For panel 38, shown 


in Fig. 18, 


0974 1s ; Z 
735 aie © = ays + 1.137 ag axial ate 


-6493(5461) + .3388(8682) + 1.137(-3785) = 2184 


or, axial panel load = 2184 1b compression. This compressive 
force, when divided by the panel and cross-section (taken 


perpendicular to the panel sweep angle ,/\p) , gives 


Oo = 2184 = ie 84 
a tw 145 x 20.09 x 273865 
Pp 
Om | 
p = 978 psi (compression) 


The shearing stress was found by taking the average 
shearing force, given for either end of the panel by the 
{at} matrix, and dividing it by the same tw used in deter- 
mining the axial stress. The shearing stress and axial 


Stress thus found were resolved by Mohr's circle to a 


37 








maximum principal stress, may, @ maximum shearing stress, 
Tmax» and the angle of the principal axis, d pe These 
stresses and angles, as listed in Table XI, refer to the 
mid-plane of the skin, since average values are used for 
their determination. These values are showm graphically in 
bere 25). 

The vertical deflections of the intermediate rib 
(station y, = 81.98) end points were calculated from the 
deflection influence coefficient matrix [C] by applying 
the theoretical vertical loads P, through Py e At the 
forward end of the rib(load point for P, ) the deflection 
was + .0601 in. At load point seven, the deflection was 


-.0647 in. These values are plotted in Fig. 26. 
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III EXPERIMENT 
1... Introduction 

In order to determine the actual stresses and compare 
them with those predicted by theory, an experimental 
investigation has been carried out on the full scale 
FSU-3 wing acquired from Chance Vought Aircraft, 
Incorporated. The experimental set-up was also designed 
to provide the Aeronautical Structures Laboratory of the 
U. S. Naval Postgraduate School with a permanent specimen 
of modern wing construction mounted such that a variety of 
test-to-theory correlation experiments could be academically 
demonstrated. 

The F8U-3 wing was structurally complete in every 
detail and was representative of current doubly tapered 
multi-cell swept wing configurations. It was divided into 
three major sections, the center section which extended 
to the wing fold rib at wing station y,, = 152.25 and an 
outer panel on either side. For this experiment the two 
outer panels were removed along with the leading edges and 
those portions of the trailing edges outside of the main 
structural box that were not integrally connected with its 
Single piece skin. Hereafter the term wing will refer to 
the main structural elements of the center section which 
consisted of seven spars, hereinafter called beams, and 


two streamwise ribs. The section was a typical thick 
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milled skin construction and the skins were tapered along 
the span as well as in the chordwise direction. 

Structural behavior was investigated with the wing 
subjected to two pure torque loads applied perpendicular 
to the elastic axis at the wing fold rib. The elastic 
axis is defined by the curve connecting the shear centers 
of all sections of the wing. Its position was taken from 
manufacturer's data as shown in Fig. 8. The torque loads 
were applied to each wing fold rib by hydraulic cylinders 
actuated from a common pressure manifold. The load 
magnitudes were measured by dynamometers connected in 
series with each cylinder and the torque was then easily 
computed knowing the fixed lever arm. Nominal torque 
values of 294000 in. lb and 336000 in. 1b were used. These 
were well within the elastic range but still of sufficient 
magnitude to give adequate strain levels throughout the 
structure for repeatably accurate measurements. 

Some electrical strain gages were already installed 
by Chance Vought Aircraft, Incorporated, at various 
locations, but for the purposes of this and future invest- 
igations additional SR-4 strain gage rosettes of the AR~-7-2 
and A-7 type were installed in particular areas. Those of 
immediate interest were located in the skins, beam webs and 


beam caps around a section perpendicular to the center 
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intermediate beam at wing station y, = 98.7. From the 
wing fold to this point the elastic axis remains nearly 
parallel to the center intermediate beam which was used 
as the reference axis of the wing. Other rosettes were 
located around a streamwise section near the root at wing 
station yy = 34.5. Here rosettes were placed on both sides 
of the very thick skin to check for differential bending 
effects. It was expected that the stresses that are 
unpredictable by simple beam theory in this area would show 
closer agreement with the predictions of the matrix force 
method of analysis. Strain gage data was programmed into 
a digital computer to obtain principal stresses and directions 
at all rosette locations. 

Deflections were measured by means of scales hung from 
the structure along the front and rear beams and read by a 
transit. Jig support deflections were also measured in this 
manner. 
2. Equipment 

The variable incidence swept wing was of a conventional 
multi-cell construction with a quarter-chord sweep back of 
42 degrees. The main structural elements are shown in 
Fig. 8. Fig. 19 is a photograph of the wing with the upper 
skin removed. Attachment was made to the fuselage at 
approximately wing station yy = 25 by two pivot lugs located 


approximately 5.3 inches aft of the main box rear beam, 
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two bumpers located at the main box front beam and two 
incidence actuators located approximately 17.8 inches 
forward of the main box front beam. The wing fold at wing 
station yy = 152.25 was accomplished by fittings with 
multiple upper and lower lugs. The upper and lower surfaces 
were single piece thick skins tapered along the span and 
chord. The front beam had numerous supporting lugs 
integrally machined into the beam giving it a very 
discontinuous variation of thickness. The structure 
primarily consisted of 7079-T6 thick skins and forgings 
with 7075-T6 sheet metal beams. 

The wing was mounted inverted on a rigid support jig 
and fastened to it at four points, the two fuselage pivot 
lugs and the two points on the main box front beam directly 
opposite the bumper points. Plywood pads were used under 
the jig to distribute the loads to the laboratory floor. 

A photograph of the mounted wing is shown in Fig. 20. 

The torque loads were applied through a fitting 
connected to the wing fold lugs as shown in Fig. 6. The 
fitting assembly was designed to provide a loading plane 
perpendicular to the plane of the wing and the elastic axis. 
It consisted primarily of mild steel parts designed to 
develop the ultimate strength of the wing fold lugs. A 
limiting torque of 420000 in. lb was selected which was well 


below the design ultimate at this station but still high 
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enough to ensure adequate strain levels in the wing. A 42 
inch lever arm was used to permit attainment of the 420000 
in. lb of torque with loads of only 10000 lb. This enabled 
the use of available 10000 1b dynamometers which were 
graduated in 100 lb increments. Instruments of larger 
capacity had 250 1b increments. 

Only the port wing was instrumented but the toraue was 
applied to each wing fold station for balance and to 
minimize warping of the structure and support jig. Each 
of the four loads were applied through a series linkage of 
attachment fittings, clevises, hydraulic cylinder and 
dynamometer as shown in Fig. 21. The four hydraulic 
cylinders were high pressure wing fold actuators and were 
simultaneously subjected to pressure from a common manifold 
to insure load equalization. Except for the dynamometers 
each linkage was designed for 30000 1b and was pre-tested 
to 15000 lb before installation. 

Load magnitude was determined and monitored by two 
separate means. On the instrumented wing two 10000 lb 
Dillon Dynamometers were used to monitor equality of 
individual loads. They were graduated in 100 lb increments 
but could be accurately read to 10 lb. The final load 
magnitudes were determined by hydraulic manifold pressure 
related to pressure-load calibration curves for the 


cylinders as found in Apvendix D. The hydraulic gage was 
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mounted at the data taking station and provided an 
excellent means of setting the load and monitoring the 
consistancy of its mégnitude throughout the run. It was 
a 1500 psi gage graduated in 10 psi increments and could 
be accurately read to 5 psi. 

The dynamometers and the hydraulic system were 
calibrated on a 300,000 lb Riehle Tensile Testing Machine 
(Appendix D). The Reihle machine had exhibited a maximum 
calibration error of only 0.35 percent over the entire 
300,000 lb range, but showed no error in the range of 
7000 1b to 9000 1b, which was of concern in this case. 
Therefore the calibration curves represent variation from 
absolute values. 

The up-load linkages were anchored to a three column 
supporting structure shown in Fig. 21. The base of this 
structure was arranged to accommodate attachment of the 
down-load linkage thereby minimizing strength requirements 
of the base elements. Considering future uses of this 
structure the base was located approximately in the plane of 
the load fitting such that single point loads could be 
applied anywhere along its reach and be readily anchored 
with a minimum of effort. 

Hydraulic pressure was supplied by an electrically 
driven Vickers V-line Piston Type Pump shown in Fig. 22. 


It was rated at 5000 psi, 1800 rpm and 1.72 cu in. per 








revolution. By means of a pressure compensator in 
conjunction with a volume limit hand control and relief 
valve, pressures could be held within 5 psi of any desired 
setting for long period of time and with very few adjust- 
ments. This feature enabled accurate repetition of the 
same loads which was essential in the data taking process. 
All strains were measured with SR-4 strain gage 
rosettes manufactured by Baldwin-Lima-Hamilton Corporation. 
In Appendix E is a list of the 485 gages, their gage factors, 
resistances and coordinate locations in the wing. The types 
of gages installed by Chance Vought Aircraft were deduced 
from resistance tests and inspection. They were all of 
the AX-5 or AR-l type. The gages installed for this 
experiment and for future multi-purpose investigations were 
either of the AR-7-2 or A-7 type. The gages of particular 
interest to this investigation are depicted in Fig. 23. 
Selected root skin rosettes are backed-up on both sides of 
the skin. Although additional back-up rosettes were 
desired throughout the inboard panels, they could not be 
installed because of time, laboratory priority and monetary 
considerations. Note that all of the gages have a gage 
factor very close to 2.0 and resistances of 120 ohms. This 
enabled the use of a single temperature compensating gage 
and factilitated reading all gages without adjusting any 


of the other instruments. 
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Selectivity in reading any combination of gages was 
made possible by routing all gage leads to a junction 
panel as shown in Fig. 24. Any number of gages could then 
be connected into either switching and balancine units or 
automatic scanning devices as future needs dictate. All 
electrical connections were soldered to minimize contact 
resistance except for the banana plug connections between 
the junction panel and the switching and balancing unit. 

A 20 channel Baldwin-Lima-Hamilton Switching and 
Balancing Unit was used. It was connected to an external 
Wheatstone bridge circuit powered by a Hewlett Packard 
Power Supply (Model 721A). A constant 6 volts was used 
throughout the tests. The bridge output was amplified, 
then fed into a voltage-to-frequency converter and displayed 
as strain in units of micro-inches per inch by an electronic 
counter. Calibration of this equipment proceeded each run 
and the method is described in Appendix F. 

The electronic counter was a Model 521DR manufactured 
by Hewlett Packard Corporation. The voltage~to-frequency 
converter was a Model DY-2210 manufactured by Dymec 
Incorporated and the amplifier was a Kintel Model 111BF. 
By using the electronic counter, readings could be taken 
faster and more accurately than with the common strain 
indicators. Strain indicators can normally be read to 


about 5 micro-inches per inch, whereas the electronic counter 
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indicates strains to one micro-inch per inch. 

Wing deflections were measured by means of scales 
hung at intervals along the front and rear beams and read 
by a transit manufactured by Keuffel and Esser Company. 
The scales were graduated in .02 inch increments but could 
be read to .005 inch. The locations of the scales are 
shown in Fig. 8. Note that a single scale could be paired 
with one of two others such that either a streamwise chord 
or a section perpendicular to the center intermediate beam 
would be defined. For instance, scale numbers & and 7 
define the test section perpendicular to center intermediate 
beam at wing station y, = 98.7 and number 8 and 11 terminate 
the streamwise rib at wing station y, = 81.98. This was 
believed to be of use in future investigations. 

Bi Experimental Procedures | 

The experimental procedures consisted of two tests, 
in each case a different magnitude of torque being applied 
at the wing fold rib. In the first test a nominal 294000 
in. lb was applied and 336000 in. lb in the second test. 
These values were initially attained by dynamometer indi- 
cations of 7000 1b and 8000 lb respectively, which will be 
the names of the two tests hereafter. The purpose of the 
7000 lb test was to establish reliability of the entire 
485 strain gages and to check for linearity of strain 


readings in the structure. The 8000 lb test was accomplished 
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to correlate the theoretical results. Therefore in the 
8000 lb test attention was focused on approximately 150 
mares of particular interest. 

The procecure in conducting the tests was straight 
forward, but because of the large number of strain measure- 
ments the tests took considerable length of time. Only 20 
gages could be read during any one run because only one 
switching and balancing unit was employed. Additional units 
were tried but the added wiring and connections, resistance 
peculiarities inherent in each unit and prolonged time 
required to adjust initial zeroes, produced an unacceptable 
drift in the zero rechecks. On the other hand load 
repeatability was very accurate and therefore a single unit 
was used and the number of runs increased. The sequence of 
operation for each run was as follows: 

1. The electronic counter was calibrated before each 
series of runs and rechecked after completion. 

2. The switching and balancing unit leads were plugged 
into 20 gage terminals at the junction panel. 

3. All strain gages and dynamometers were zeroed. 

4. The hydraulic loads were applied and adjusted by 
reference to previously known values of hydraulic pressure 
and strain gage readings from the same check gage used in 
each test run. 


5. The dynamometers were read and compared with each 
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other to insure that equal couple loads were being applied 
to each wing. 

6. The strains were read and recorded twice. 

7. The loads were removed. 

8. All zero values of strain were re-checked and 
recorded. 

Extreme care was employed to achieve load and strain 
repeatability. The electronic counter was calibrated 
before each run and checked afterwards. A maximum departure 
from linearity over a range of 1000 micro-inches per inch 
was only 5 micro-inches per inch. The average departure . 
was closer to 2 micro inches per inch. This effect on 
maximum strain levels encountered (approximately 200 micro- 
inches per inch) was therefore quite neglirible. 

In addition particular care was exercised to minimize 
7 strain gage zero drift. For instance, the variable 
temperature effects of sunlight and outside electrical 
interference with the sensitive instrumentation were 
virtually eliminated by conducting all tests at night. At 
least one hour warm-up time was given all instruments to 
insure steady state conditions. All initial zeroes were 
re-checked three times before loading. All permanent 
electrical connections were soldered to minimize contact 


resistance and, as a double precaution, a test was repeated 
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if any electrical leads had been disturbed or if the 
voltage supply varied any time during the run. By reading 
only 20 gages per run the time per run was minimized, 
thereby reducing the critical time the system would be 
exposed to outside effects. Two strain readings at load 
were taken also as a double check. With such care, average 
zero rechecks within 1.5 micro-inch per inch were 
continually demonstrated. 

Two strain gages were selected to act as cross-checks 
on the repeatability of the loads, results and test 
conditions. Gages 32 and 53 were monitored in each test 
run. The magnitudes of strain on each runwere then 
compared to the original values. In over 40 runs the 
average departure from mean values of the check gages was 
only 1.39 micro-inch per inch. The mean value for gage 
53 in the 8000 lb test was 178.6 giving an error of only 
0.78 percent. Since the hydraulic pressure reading 
required to produce this strain level was always the same 
the load repeatability was also within an accuracy of one 
percent. 

Accuracy of measuring the absolute load magnitude was 
then depencent only on the accuracy of the dynamometers and 
hydraulic pressure gage. During the weeks of testing, 
difficulty was encountered in re-setting the dynamometers 


to the same zero. Also there was noticable differences 


50 








in readings under load conditions which verified suspicions 
of unreliable zero settings. On the other hand the 
hydraulic pressure gage readings checked consistantly with 
strain sage 32 and 53 readings and therefore the pressure 
gage was felt to be more reliable in establishing magnitudes 
of applied loads. 

Deflection measurements were made for only the 8000 
lb test condition. Since this loading condition and 
consequently the deflections were repeated over 20 times 
it was deemed necessary to only sample the deflections 
periodically throughout the test series. Although more 
deflection measurements were taken to ensure repeatability 
in the series only three complete sets of data were 


recorded. 


4. Experimental Results 

From the recorded strain data the only calculations 
required were those involved in determining the principal 
stresses, maximum shear stress and principal axes at each 
rosette location. These calculations for all rosettes in 
the 7000 lb test and those particular rosettes of interest 
in the §000 lb test were performed on a Control Data 
Corporation 1604 Digital Computer. The results along with 
the FORTRAN program and the measured strains are found 


in Appendix G. 
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The cata for the 7000 1b test was pertinent to this 
experiment only to show linearity of strain readings. A 
sampling of net strain readings from various locations for 
both loading conditions is shown in Table XIl. Satisfactory 
linearity is demonstrated by noting that the average 
difference between measured strains for the 8000 lb test 
and those extrapolated from the 7000 lb test is 2.3 micro- 
inches per inch regardless of the magnitude of the 
measurement. The discrepancy is obviously quite small but 
it would be still further diminished if the accuracy of 
the applied loads were considered. Therefore it can be 
concluded that linearity of strain readings was indeed 
achieved. 

Another result of the 7000 lb test was the discovery 
of 9 faulty rosettes. However none of the faulty rosettes 
hampered the results of this experiment. The faulty 
gages are indicated in Appendix E. 

Calibration of the hydraulic system after the completion 
of the tests revealed that instead of 8000 lb, actually 
8400 lb loads were being delivered to the linkages. 
Considering a gage reading accuracy of 5 psi and the slope 
of the calibration curve as 10.84 1b load per psi hydraulic 
pressure, the accuracy of the applied loads would be plus 
or minus 50 lb or 0.6 percent. The difference in desired 


and achieved magnitudes is not of real importance however 
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since linearity of strains enables interpolation to the 
desired values. 

Since a degree of uncertainty always exists concerning . 
the extent to which a pure torque loading is realized, an 
equilibrium check was made at the instrumented section 
perpendicular to the CIB at y, = 98.7. At this section 
all rosettes were oriented such that the diagonal gages 
were at a 45 degree angle to the section, which facilitated 
shear flow calculations. Therefore theoretically 
calculated shear flows could be readily checked with 
experimental results since the theoretical calculations were 
all ready completed in the determination of the loading 
for the idealized structure as shown in Appendix A. 

The experimental shear flows were calculated by 
considering only the strains in the diagonal gages. This 
essentially eliminated the effects of any normal stresses 
caused by unknown bending loads and the results would 
better represent the shear flow induced by only the 
effectively applied torque. From the experimental data 
for the 8000 lb load, shear flows were calculated using 


the relation, 


qa= Etee (15) 
1 + Y 


ws * 
° ee 
agdo* 
= 
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These results are listed in Table XIII and are labled on 
the schematic diagram in Fig. 25. The values of shear 
flow in the top and bottom skin of the same cell compared 
closely. The average value of each cell was used as also 
shown in Table XIII and Fig. 25 and compared to the 
theoretical values. The theoretical shear flows for 

T = 336000 in. 1b were 3.85 percent higher than the experi- 
mental values. The total torque resisted was then computed 
as the sum of the (2Aq) for each cell. This value was 
323,785. in. lb which is 8.2 percent lower than the 
352,800 in. lb possible with an 8400 1b force acting on a 
42 inch lever arm, and 3.6 percent lower than that used in 
the analysis. 

Neglecting the accuracy of the experimental loads and 
shear flow determination this means that 720 lb of the 
applied load of 8400 1b was not effective in producing pure 
torque. In this regard, during the tests there had been 
visible evidence of streamwise twisting of the wing fold 
rib and a measured spanwise tilt of the entire wing. There 
was also an element of uncertainty regarding the exact 
angle of the elastic axis since it was scaled from a 
reproduced drawing. Coupled with understandable inaccuracies 
of maintaining a constant loading plane, perfectly parallel 
lines of force and a perfectly rigid loading jig, the amount 


of the bending load component is within reason. Although 


54 








no calculations were made to determine the magnitude of 
the bending component its existance can be qualitatively 
confirmed by observing the magnitudes of the strain 
readings for the perpendicular gages at the section being 
analyzed, namely rosette numbers 131, 46, 49, 52, 55, 58, 
etc. 

Using a tonque,valuewef 3234755,.an. 1b the theorewncel 
shear flows were obtained by interpolating between values 
computed for the 294000 and 336000 ins Jb condition. 

They are also shown in Table XIII and Fig. 25. As would 
be expected the comparison with experimental values is 
very good. The average error in the theoretical shear 
flows is only 1.95 percent. Considering the loss of 
effective applied toraue the torsional equilibrium check 
was believed to be quite close. 

The computer program described in Appendix G solves 
the well known rectangular rosette equations for d,.,, 
Pes 7 was and Dp. The input to the program was the 
value of strain for each gage in the rosette. The lowest 
numbered gage was consistantly called €,, the diagonal 
and perpendicular gages were the next higher numbered gages 
respectively. The principal axis was then computed with 
respect to the €}, axis. To compare these results with 
theoretical values the orientation of the principal axes 


were adjusted to the yw axis of the wing by simple arithmetic. 
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The root area rosettes used to compare with results of 
the theoretical root analysis were numbers 13, 19, 28, 367, 
279 and 161 on the inside surfaces and back-up rosettes 
numbers 481, 478, 472, 466, 457, and 484 on the exterior 
surfaces. Differential bending was readily apparent in all 
of the root panels as seen in Table XIV by noting the 
differences of strain in matching gages on opposing surfaces. 
The effect is less in panel number 21 which lies near the 
elastic axis. The average of the two gage readings was 
taken as representative of the strain at mid-thickness 
of the skin. The values of Omax? eae and Dp were then 
hand calculated using the equations in Appendix G and are 
listed in Tables XI and XIV and plotted in Fig. 26. 

The deflection measurements were of concern primarily in 
establishing test conditions. In this regard the average 
values of the tota. deflections are plotted in Fig. 26. 

i? ais clearly shown that the support jig rocks forward and 
the entire structure tilts right-wing-down. It was also 
noted during the test that the amount increased slightly 
with each cycle which indicated plastic yielding of the 
plywood pads under the support jig. Obviously any tilt 
would adversly affect the lines of action of the applied 
loads since the load linkages were positionec to accommodate 


only one loading plane. This would account for some of the 


discrepancy in achieving pure toraue loads. It would be 
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wise to consider adjusting the load fittings and linkages 
such that a more symmetric loading could be realized in 
future tests. Furthermore it would be advantageous to 


replace the plywood pads with a more suitable material. 
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IV COMPARISON OF RESULTS 
Extreme care was taken to ascertain strain and load 
jwze.vuucs and insure repeatability of loads to minimize 
any experimental error in the determination of the compar- 
itively low stress levels found in the critical root area. 
All of these errors were kept well below one percent and 
logically do not appreciably influence any comparison with 
theory. 

The only experimental error that would noticeably 
affect the final results would be the understandable 
discrepancy in developing pure torque loading on which 
the theoretical analysis was based. The equilibrium check 
at the section perpendicular to CIB, at y, = 98.7 indicated 
that the torque developed was 3.6 percent lower than that 
used in the analysis and nine percent lower than actually 
applied at the tip. This meant that about 8.6 percent of 
the applied load introduced unwanted bending effects that 
were not accounted for in the analysis. The cause of the 
torgue discrepancy is correctly attributed mainly to the 
method of loading and tilting of the support jig. Other 
than reducing the accuracy of the experimental results 
to a comparatively minor degree the effect of the small 
bending loads would not measureably influence the 


comparison with theory in such a complex structure. These 


factors become even less important in establishing the 
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general validity of this analytical method to an actual 
complicated wing when one is reminded that the literature 
mentions 100 percent disagreements in comparing results of 
other methods applied to simplified thin-skinned laboratory 
specimen. 

Experimental and theoretical results for the root area 
are listed in Table XI, Three observations may be made 
immediately. The sense of maximum normal stress, the 
orientation of the principal axis, and magnitude of stresses 
at the critical rear panel area agree rather well. Maximum 
shear stresses differ by only 25 percent, and although one 
of the maximum normal stresses differs by 40 percent the 
other is only seven percent. On the other hand, the 
theoretical stresses obtained at the leading edge are much 
too high. This is shown graphically in Fig. 27. Since the 
equilibrium checks made during formulation of the problem 
indicated excellent equilibrivu.: existed, and agreement of 
theoretical and experimental results at the rear of the 
wing is good, the leading edge divergence was likely to be 
the result of an initial assumption. 

Attention was immediately directed to the root boundary 
conditions. It was assumed that the pivot rib (considered 
to be the root of the idealized model) was infinitely stiff 
and completely resisted chordwise deformation. This is 


still felt to be a reasonably accurate assumption. However, 
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re-examination of the drawings of the pivot rib indicates 
that the rib width varies by descreasing from rear to front 
substantially. Moreover, the skin thickness varies the 

same way. It would therefore appear that significant twist 
about the streamwise axis occurs towards the front of the 
pivot rib. This has the effect of removing the root 
restraint imposed upon the theoretical analysis, and allows 
stress relief to occur towards the leading edge in the actual 
case, 

The sizeable differences in back-up gage readings 
vould further stimulate interest into investigating the 
plate bending and twisting effects of individual panels 
about their respective axes. Unfortunately, insufficient 
back-up gages were installed to establish any definite 
conclusions, but examination of the stresses at the 
available back-up locations and the stresses measured by 
neishboring rosettes on only one side of the skin evokes 
considerable concern in this area. Panels 25 and 30 
would be particularly well suited for such an investigation 
in the future. 

It is therefore believed that the analysis method 
can achieve accurate results for this wing. The boundary 
condition in question could be removed entirely from the 
analysis by moving the theoretical root to the aircraft 


center line. This would introduce only 16 new unknowns 
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to the analysis, all within the additional rib required 
at station hoe 25.095. The addition of 16 unknowns would 
in turn require conversion of the computer program to 
machine language in order to accommodate the increased 
storage requirements. This could be accomplished by use 
of the Fortran MAP program available at the U. S. Naval 
Postgraduate School Computer Facility. 

There were no load attachment points on the specimen 
corresponding to the analytical load attachment points. 
This could have enabled direct experimental verification 
of a greater portion of the flexibility influence coefficient 
matrix (C| (Table VIII). However two node points can be 
compared here. The measured net vertical deflections at 
scale numbers 8 and 11 (Fig. 26) are compared to deflections 
caused by the analytically calculated vertical loads P, 


through P At scale number 8 the theoretical value differs 


from the ocean value by only 16.8 percent and at scale 
number 11 the difference is only 14.8 percent. The better 
agreement between theoretical and measured deflections was 
expected(Ref. 9). 

Improved accuracy in the solution could undoubtedly 
be obtained by reducing the size of the structural grid. 
This would at the least double the number of unknowns. 


With the additional element flexibilities required, 


computer programming should be extended to do more of the 


61 








labor. As has been pointed out by Rattinger and Gallagher 
(Ref. 2), up to three man-years may be reauired for 
completely progremming a solution using a displacement 
method. It is felt that the force method of this report 
could be extended to-twice the current size with about one 
man-month additional effort, if project familiarity were 


equal to that of the writers. 
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V CONCLUSIONS AND RECOMMENDATIONS 
a Conclusions 


It can therefore be concluded that a valid comparison 
was achieved between experimental and theoretical stresses 
as predicted by the matrix force method of analysis, using 
a structural idealization proposed by Wehle and Lansing. 
Further, this analytical method gave remarkably reasonable 
results in the critical rear root area of the highly 
complex structure. 

The failure of the method to agree with actual stress 
conditions in the forward root area is attributed primarily 
to the simplifying assumption that the pivot rib was 
infinitely rigid. 

Le Recommendations 

It is recommended that all cover skin rosettes, 
especially those between the intermediate and pivot rib, be 
backed-up. This would greatly enhance the experimental 
potentiality of the laboratory by providing means to 
investigate the plate phenomena of bending and twisting of 
individual panels about their respective axes. 

In addition greater academic value could be realized 
from the analytical results if means were provided for 
applying single point loads at the intersections of the 
intermediate rib and the beams. This would permit further 


interesting academic demonstrations by utilizing the matrix 
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of flexibility influence coefficients produced in this 
analysis. 
It is further recommended, to improve the accuracy 
of the analytical predictions along the entire root area, 
that the root boundary condition assumed in this analysis 
be eliminated by moving the theoretical root to the wing 
center line rib. This would not require substantial 
increase in analytical complexity and would greatly 
enhance the agreement between theory and actual measurement. 
It is finally recommended that the experimental 
method of loading be improved to enable better development 
of pure torque. In this regard an exterior rib, clamped 
around the section perpendicular to the center intermediate 
rib at Wing station y, = 143, would eliminate the stream- 
wise twist of the wing fold rib. Regardless of the method 
of loading, the plywood pads under the wing support jig 


should be replaced with a more suitable material. 
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TABLE IT 
COMPONENTS FOR APPLIED LOADS 







ds ani 
Cell Shear Flow 
ge stn" 











3581. 83' 4781.46 


3790.74 
4815.86 


265.74 99. 84 


















4869.70 
5920.33 
6330.00 


IS ao 
149.82 


302.03 
398.78 
457.81' 






















5335.69 
5310.11 


171.99 
a itieate 
166.23 








6128.50 








“472.81 
442.47 

















4840.75 5476.53 





*e= (90-2/\) 
(See Appendix C for resultant {P} matrix) Ee 
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1 
(Ap sinh,) 


3223.13 
3056.91 
3443.56 
3405.79 
3059.59 
2561.21 
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TABLE III 
THE EDGE SHEAR FORCES ON THE PANELS 








A 
Wy “y ao q 
y - io 
age 
BP: 
1g 
ela ie Te eee ee 
in. in. ae - % la/w 
1 14340883 |14.328741 | 2.102769 | .826979 | 2.100988 
2 , . 13.504185 | 1.815557 | .916666 1.815079 
B - : 13030300 1.701104 .971279 1.701084 
4 : 12.575467 | 1.667853 }1.015915 1.667834 
5 12.137011 |12.129637 | 1.719123 |1.067989 1.718079 
6 11.804608 | 11.798490 1.120635 
if 79.3021 «779006 
8 74.9206 | .762864 
9 71.3846 770186 
10 68.5890 . 789691 
Al 66.4567 «819565 
12 . 848558 
13 (f ; . 875000 
15 20. 870584 179.3021 . 686557 
17 19.099132 |74.9206 «707058 ' 
19 17.787945 {71.3846 «732535 
al 16.511403 168.5890 .761623° 
| Be 15 .239329 166.4567 ~ 795943 
25 14 .141857 |64.9038 . 8342905 
28 14 .174828 164.9038 6832787 
30 15.245680 166.4567 ~ 795767 
32 16.513296 168.5890 6761544 4 O244y4 
34 17. 787956 |71.3846 ° 4.013086 | .732561 3.85592 | 
36 19.102224 174.9206: 3.922088 | .707129 3.73698 
38 20. 887886 |79.3021 3.796559 | .686564 3.58680 
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-.019508 
~.012482 
-.175678 
~. 136654 
~.098466 
-.059240 
-.017702 
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oe 
-000001 


-000001 
- 000833 
-000650 
-000558 
-000511 
s00052 | 
-000584 
-000004 
~000004 
- 000004 
-000003 
-000003 
-000003 
~000003 
-000003 
- 000003 
-000003 
-000004 
- 000004 
-000004 
~000004 
1.000024 
-.005421 
-.003903 
-.002792 
-.001859 
-.000979 
-000979 
~001859 
-002792 
- 003902 





0667642 
- 667642 
- .275926 
188647 
142422 
-115168 
110933 
» I2eaG 
1.926192 
~ 375191 
o 349377 
2 324341 
e299KK4 
0278411 
277603 
299146 
0324117 
0349175 
374914 
1.924836 
000000 
000000 
000000 
093957 
0216014 
2545211 
983284 
486114 
~485306 
~982178 
-1.543881 
-2 «214480 
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481468 
~.481468 

554455 
~. 368351 
~.279407 
-.227115 
-.219338 
-.245798 
1.389068 
13270320 
- .658480 
-.611294 
- .564370 
-.524728 
- .523376 
~ 563992 
-.611070 
-.658313 
1.270295 
1.388115 

000000 

000000 

~000000 
1.977363 
4.253180 
2.967102 
1.889707 


— .935407 


-.934056 


-1.887978 
-2 965149 
-4.251059) 
-1.975265. 
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me. 36 37 38 oe 40 oa 42 
2338325. .224309  .115240 .002875 003653 + .003226  .00215h 
meee so5) -.224309 -.115240 -.002875 ‘-.003653 -.003226 -sooeTSE 
2420912 .276236 .139517 .-.006221 -.007048 -.006581 -.0043C% 
2469520 .306469 1154475 .-.007239 -.009888 -.008097 -.005116 
~.424733 2374470 .188657 -.007936 -.011660 =-.011414 -.007127 
~.345168 -2498116 =. 252915 -.008626 -.013178 -.013967 -.010959 
Mssi2 ~-.482222 -.662N07 -.009387 -.014524u -.015766 -.013093 
~.373523  -2540398 -.742320 -.010519' -.016277 2017668 -.014673 
976090 0 6471K5 0332476 - 008294 -010538 - 009308 ~006215 
6892646 591822 .304054 ..007585 .009637 .008512 .005684 
4831230 .55110% .283124 .007063 .008974  .007926 .005293 
~.930730 6511612» «262845 = 6006557 =. 008331 =. 007358 = «004914 
~.859285 -1.184810 .242668 .006054 .007692 .006793 .004536 
~.798928: -1.101588 -1.478073 .005623  .007151 .006316 .004218 
~.797021:-1.099057 -1.474699 .043634 .083555 .121239 .156745 
~.858872 -1.184347 .242500 .047020 .090080 .130648 .168909 
~.930564, .511519 .262742 .050945 .097555 .1K1554 -.065446 
6831290 .551066 .283055 .054883 .105098 -.106843 -.070505 
6892568 .591687 .303920 .058929 -.1584708 -.118719 -.0757C2 
2975353, 646566 =o. 332108 «=. 214175 -.169057 -.125359 -.082724 
000000 .000000 .000000 -.845400 -.682561 -.511949 -.338279 
.000000' .000000 . .0000CO0 -.101080 -.151244 -.160424 -.132661 
‘-000000 ~000000 ~O000CO -—.272641 -.496083 -.679040 mae 
1.228755 .743506 .360409 .036391 .057173 .061430  -.047151 
2.884878. 1.836419 .917554 .032660 .053981 .057958 .044993 
4.514167 2.908448 1.463261 027402 .046122 052098 041572 | 
2.874906 4.044536 2.040715 .020722 .035005 .040418 .0345S8! 
1.423008 2.004500 2.718373 .011957 .020057 .023218 .020309 
-1.421101 -2.001970 -2.7150C0 .026037 056322 .091668 132169, 
=2.872586) —4%.041543 -2.037511 .058223 .123693 .198278 .2821921 
-4.511682, -2.905548 -1.460161  .095934 .201806 .320804  .204866 
~2.882332) -1.833556 -.914522 .138535 .290150 .200233 .125686 
-1.226285, -.740776 -.357518 186278 122717 073609 ~° .042194 | 
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43 Ls Ls; Le Lg 8 Lg 
5001012 -.054604 —.043469 -.033989 -.024259 -.016002 -.007931 
~.001012 .054604 .043469 -033989 + .024259 .016002 #£.007931 
~.002289 -.060804% -.057410 -.048118 -.038192 -.025181 -.012564 
-.002595 ~-.0544H17 -.053040 -.046621 -.036420 -.023823 ~.011803 
~.003516 -—.051717 —.051458 -.046558 -.038353 -.024920 ~.012277 
~.005332 -.051323 -.051813 -.047802 -.040717 -.028631 -.014035 
~.007721 -.054078 -.054891 -.051005 -.043961 -.031721 -.016898 
-.008653 ~.060603 -.061513 =-.057158 -.049264 -.035548 -.018936 
002921 2105861 -e125411 -.098059 -.069989 -.046168 -.022880 
2002671 .096811 .140276 -.089677 -.0640C6 -.042221 -.020924 
002487 8.090150 = .130625 2165678 -.059602 -.039316 -.019485 
.002309 .083690 .121264 .153806 .185426 -.036499 -.018089 
.002132 077266 .111956 2141999 6171192, 6195711 -.016700 
.001982 .071839 .104%092 ~132025 .159168 .18196%  .202870 
2190202 6071677 =. 065837) =— 6.055302 = .0438702=—S 029019 «w. 04 187 
-.029836 .077240 .0709%6 .059593  .O47274 .031271 .015288 
-.032326 .083687 .076868  .064568 .051220 .033881 016564 
-.034826 .09015% #.082811 069560 -055180 #£.036501 2017845 
~.037393 .096803 .088916 .074687 «= .059248 .039192 .019160 
-.040861 .105781 .097162 .081614 .064743 .042827 .020937 
=-.166714 1.000056 .845400 .682561 .511949 .338279 .16671% 
-.076154 -.809454 -.956900 ~-1.048542 -1.0934C07 -1.094728 -1.060437 
-.925501 -.809454 -.785339 =-.703703 -.574792 -.406460 -.211091 
2026497 .059017 -.136139 =.081876 -.0445S4 -.023150 -.008466 
-024997 045401 -—.1000TS -.258866 -.177913 -.111075 °-.052196 
' 6023065 .042970 -.059596 +.172438 -.299429 -.190891 -.091747)/: 
2019508 .040403 -.024154 -.096280 -.177961 -.268948 -.130310 
2012482 .029481 -.001239 -.036442 -.076770 -.122486 -.171152 
.175678 -.029643 -.037C0%. -.040256 -.038491 -.030409 -.017471 
2136654 -.040591 -.055083 -.062795: -.061173 -.048327 -.026295 
2098466 -.043161 -.064CKH1 -.075881 --.073921 -.056011 -.030209 
2059240 -.045587 -.071475 -.085651 -.080715 -.060050 -.032451 
2017702 -.059212 -.086841 -.090858 -.066614 
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a 


16 
18 
20 
Be 
Qu 
ey 
oil 
BS 
5D 
37 
Be 


2nd post 
3rd post 
4th post 
5th post 
6th post 


7th post 
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TABLE VII 


CHECK OF CJ EQUILIBRIUM 


3.58977(a5,) 
3.73758(ay9) 
3.85593(ajg) 
4.02490(aj7) 
4 .25897(ay6) 
4.578887 (aj 1) 
4.359051 (ay 3) 
4,153562(az>) 
4 .013086(a, 4) 
3.922088(a, ,) 
3.796559( ag) 


.826979(—.3274394) + 1.0 = .8091530 
-916666( -.1468006) 


3.922723(aj9) + 
4.013088(a, 9) + 


4.154038( a) 7) + 
4 ,360868( a6) + 


4.024WW(ay>) = 


3.85592( az) - 
3.73698(aj,) = 
3.58680( ag) = 
+ 10.95333(ao) = 


©971279(-.1547262) + 1 


1.015915( =. 142922) 


+ } 


1.067989(-.1018803) + 1 


1.120635(-.0092286) + 1 
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- «7766476 
- «7840807 
- .8544135 
-1.0103419 


Beg 
AU 
=a 
aya 
AY 
Ang 
2H 8 
a7 
AYE 
YS 
ayy 


- 0 


- 0 


+ 


+ 


.00000018 


* +,00000028 
+.00000033 


‘=.00000077 


00000053 


- +,00000082 
+, 00000047 
~=.00000068 


‘=.00000026 


-000004.30 


- .0799385 
+.0795013 
+.0730701 
+.0707227 
+,0367795 
+ 0.00000 
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-859 .SU667TECY 
~3524 62157683 
~3361.354572430 
~3227 .87656528 
~2977 .4u3 32845 
~2714.088 87257 
767279220127 
904.957 24172 
592.9335C626 
~17.28751€87 
~720.49754617 
€552.00505352 
8725 85882282 
78702685 E331C 
4969 665661299 
1010.45895720 
5460.770 78527 
8681.69078517 
8356.9425C321 
7980.7227118C 
4504214537561 
769.2265845€ 








: TABLE XI 
| RESULTS AT STATION yy = 34.5 


Comparison of Maximum Normal Stress, Maximum Shear Stress 
| And Principal Axis Deviation From Y-Axis. Left Wing, 
| Forward Of Y-Axis Is Positive. Theoretical Loading 
| Corresponds To Experimental Loading Of 
| 336000 in. 1b Torque 








Maximum Shear Stress Principal 
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retical mental retical 
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TABLE XII 
COMPARISON OF LINEARITY OF STRAIN READINGS 


AT VARIOUS LOCATIONS FOR THE 7000 LB AND 8000 LB TESTS 


ae ; Strain ( V in/in) 
nera cation age a 
Fe iirie | | 8000 lb Load 














Number ! 
Extra— Differ— 
8000 Lb polated ence 
test | from 7000! Between 
lb test j{ Measured | 
And Extrar 
Upper Skin : | 
Outboard Near es: =~ 5e - 52 - 59 7 
Leading Edge DBS +123 +147 > +141 6 
ILS + 31] + 4y » + 35 9 
Outboard Near’ 43 f = 21 - 22 ; = 23 1 
Trailing Edge yy i + U6 + 53 f + 53 0 
45 | + 2 + 4 : + 2 2 | 
Inboard Near 16 - 11 - 12 - 13 od 
leading Edge lsh | mmc + 83 | + 82 1 | 
18 eee + 30 [ONE Be 2 
Inboard Near 89 + 8 + 8 1 + 9 1 
Trailing Edge 90 | + 50 fos | +57 L 
91 } = at = f + 13 y 
Lower Skin | 
Outboard Near 169 | +755 | + 66 + Gs 3 
Leading Edge - 170 +141: +162 +161 . 
ify ~ 36 - 41 ~ 41 0 
Outboard Near 394 eye e320 ="30 0 
Training Edge 395 +150 +170 +171 1 
396 n i + 49 ; + UG 3 
Inboard Near 161 ae + 17 + 14 3 
Leading Edge 162 i +10 + 15 f+ 11 4 
i 
Inboard Near 367 | aS + 14 +15 i 
Trailing Edge 368 - 42 - 50 + 48 2 
- 369 + 21 + 2y + 24 0 
FB Web 152 - 7 - 10 ae 2 
— 153 alloy +176 Bot EYAS, 3 
154 ee) + 25 + 29 4 
CIB Web 273 + 2 + 3 + 2 1 
274 =~ 25 © =- 28 - 29 1 
: 275 - 2 - 1 - 2 ils 
RB Web 4Q9 ~ 20 - 22 - 23 1, 
410 +134 a Sole gol 8, 2 
a sand | + 5 + 8 3 
Average Difference Between. Measured and Extrapolated = 2.3 


Q4 








Celliposette Twice 


No. 






Web 152 
a 131 
; 169 
Web} 181 
196 

b 46 
216 

Web 225 
25/ 

Cc 4g 
252 

Web 273 
261 

d De 
291 

Web; 300 
312 

e 5 
Se 
Web 342 
Spy 
4 58 
394 

Web 4Q9 


86.2 


105.8 


27.6 


140.0 


143.4 


138.0 






TABLE XIIT 


EXPERIMENTAL AND THEORETICAL SHEAR FLOWS AND TOTAL 


TORQUE REACTED AT THE SECTION PERPENDICULAR TO 





CIB AT Yw = 98.7 




















| Shear Flow,a,(Ib/in.) ss CTS Torque 
| Fxnerimental | Theoretical | Exp. | 
Initial | Average | 323,785 | 336,000 
8000 1b | 8000 1b] in. lb vB vamealle test 
test load 
222% eral 206.2 213.95 
2d 211 206.2 213.95 | 18188.2 
2a) 
224 200 232.5 241.35 
5 
415 415 438.7 455.30 | 43907.0 
445 
48,8 51.3 49.9 Bie o 
Denice 
495 488 488.6 507.08 | 62268.8 
481 
| 
10.9 10.9 4.7 4, 83 
Bad Gage 
Bie - 509 493.3 71260.0 
505 
5.8 p Veo 
(pene 
543% 478 475.5 68545 .2 
478 
sina5 38.4 40.7 
Mi. 3 
414 432 434.8 59616.0 
450 
45 * 430 434.8 451.19 


-* Considered Stray Values 
f 
| 

{ ‘ 


03 


= =e 


17774.4 


46239.0 


62345.4 


69062.0 


68186.7 


; 60002.4 





1842.1 
48170..3 
64703.9 
71667 .4 
! 
| 


70752.4 } 


| 
| 


| 
62263. 9" 


! 
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EXPERIMENTAL RESULTS AT y,, = 34.5 





Tmax |  Dytt 
psi 
| 759 nseyAce tit : 
690 +4906! 
| O 
468 +14, -46, -~16 -576 572 -4yr sur 
467 -49, -95, -/2 
466 +24, +7, +15.5 
u7y 22, +28, +3 
473 +136, +78, +107 +868 800 —44o35' 
472 +25, -13, +6 , 
480 -6. ~16, +11 +555 525 ~ 39°26! 
| 479 +72. +64, +68 
478 +26, +4 +5 
483 -78, -66. -12 -800 535 ‘| +18939! 
482 +59 .5 +59. +59 
4u8l +35, +39 , io, | 





* Two gage rosette. 


** Angle measurec from sweep angle of CIB. Plus angles are measured in the 
direction of the diagonal gage of the interior rosette. 








F8U-3 Wing With Upper Skin Removed 
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Fig. 20 


F8U-3 Wing Mounted Inverted 
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Vickers V-Line Piston Type Pump Assembly 
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Fig, 24 
Instrumentation 


Hydraulic Pressure Gage 

Strain Gage Lead Junction Panel 

20 Channel Switching and Balancing Unit 
Wheatstone Bridge Circuit 

Power Supply 

Electronic Counter Assembly 
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Pig, 2D 


EXPERIMENTAL AND THEORETICAL SHEAR FLOWS 
AT CUT SECTION PERPENDICULAR TO C.1.B. AT yw = 98.7 
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BXPRRIMENTAL RESULTS: 
SOOOMIb Test (T = 324/765 dns 1b) 


450 S16, «S05. 480 415 » 21 
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Direct Results From Rosette Readings. 
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Final Results After Averaging Skin Values. 
_ THEORETICAL RESULTS 
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FRONT AND REAR BEAM DEFLSCTIONS 
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APPENDIX A 
THEORLTICAL DETERMINATION OF SHEAR FLOWS 
AT SECTIONS PERPENDICULAR TO CIB AT 
as 98.7 AND Yi yi 

The theoretical shear flows at sections perpendicular 
to the CIB at y,, = 98.7 and Y= 74.3 were calculated for 
two reasons. First the shear flows at both sections would 
be used to establish the theoretical loads at streamwise 
rib, y,, = 8498. Secondly, the shear flows at the section 
perpendicular to the CIB at y,, = 908.7 would be determined 
experimentally and comparison with theoretical values 
would facilitate a check of the degree to which pure 
torque loadirg was achieved. Both cross-sections are 
shown schematically in Fig. Al. The thicknesses were 
average mid=-panel values and the web heights were taken 
between the mid=panel of the upper and lower skins. 

A box structure with several cells wil] have one less 
redundant than the number of cells. In this case there is 
five redundant webs. It is desired to write six equations 
in the six unknown shear flows. This was done by equating 
the angle of twist of one cell with the remaining five, 
which gives five equations, and then writing an equilibrium 
of torsional moments equation. 

The first five equaticns were obtained by equating 


the angle of twist per unit length, @, of one cell with 


94 








Pee a ef en ee ee en ED | eee 


the other five using the well known expression for a 


box beam, 


gis | 
—_ \ “rg =e ee (A4 j 
[2A tG 


The equilibrium of torsional moments may be written 


as 


= 
T=) ee ee (A2) 


In using equation (A2) the summation is carried out 
around the entire perimeter of each cell. The unit length 
L and the constant factor 2G drop cut leaving, for cells 


a and b, 


q. As - » qs (A3) 
Lika t es 
The value of q for any exterior web of cell a is q,; and 
for the interior web is (qa+ ~ Qp4). Using the 


abbreviations, 


Oaa = ) & 

ce ib 

Ge ) As 

iE 

jie As] 
ae _ t] a--b 


oD 





equation A3 is rewritten. 


& ae C 
dat Oaa - Wt Dab = Abt Opp %at Car (ayy 


fe A 


a Ay Ap 
The terms On. andOy, represent Summation around the entire 
perimeter of their respective cells and Oo» the value for 
the interior web. 

For each cell the values of O were easily calculated 
and are shown in Table Al. The enclosed areas were taken 
as the average web height times the distance between webs. 


Bquating the angle of twist per unit length of cell 


a to the remaining cells gives the foilowing equations: 


| c 
i dat Oaa 7 At One 

















a | 
ie q S 
ape. | Apt Opp = Gat Dap = Ast One 
An 
= Get One = dpe Une = Ade Sea] 
C 
| 5 (A5) 
aa | dat Oaa ~ Get Yea ~ Get Vde 
d 
i . = . 
= a | ee om e q+ CAE WPt O oe 
a | 
oa 7 iS 
So ae Ove tet i 
in 


Substituting values of O from Table Al, performing the 
indicated arithmetic and rearranging, one gets five equations 


in terms of the six shear flows with constant coefficients. 





These are shown in Table A2 for both sections. The sixth 
equation is also shown in Table A2, in which case the 
coefficients are simply two times the cell areas. 

The solution to these six simultaneous equations was 
obtained on the CDC 1604 Digital Computer using a FORTRAN 
program with a CO-OP identification of F2 UTEX LINE QN 
march is anciuded herein as paresmrO2 thmgush 11Gme It 
utilizes Gauss’s method of elimination with row pivoting 
and back substitution and is designed to give solutions 
for one or more column vectors forming the right side of 
the set of equations. This means it would pive solutions 
for cne or more values of applied toraue in this case. 
Three torque values were used as shown in Table A3. The 
vaiues of 294000 in. lb and 336000 in. 1b were used in the 
analysis. The value of 177825 in. lb was used to check 
the computer results, The check values of gq were hand 
calculated using a iterative procedure sugrested by Bruhn 
(Ref, 11) and within slide rule accuracy agree well with 


the computer solutions. 


a 
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SCHEMATIC DRAWING OF SECTIONS 
PERPENDICULAR TO THB C.1.5b. 








Cii.B. F.B. 


Yy = 743 
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TABLE “Al 
SECTION PROPERTIES 
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Section y,, = 98.7 | Section yee 74.3 
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178.250 







250.743 369.016 
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119.206 












327 . 360 250.572 


78.556 





211.850 162.757 
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TABLE AIII 
SHEAR FLOW RESULTS FOR VARIOUS TORQUES 
AT SECTIONS PERPENDICULAR TO CIB AT 


Vw = 98.7 AND Yuaau 74.3 







Shear Flows Weyagayr 
‘Hand 


Calculated Computer Results 



































in. in. "Ib "Ib tay 
134.28 135.289 223.675 255.628 | 
| a66.09 | 167.416 | 276.790 | 316.332 | 
196.21 | 197.292 326.185 372.782 
222.05 222.450 367.779 420.319 
229.65 229.231 378.990 433.132 
203.48 201.905 333.811 381.500 | 
| 





opee 02 
398. 384 
443.699 
4Y7 921 


213.945 
455.296 
507.084 
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451.188 


431.719 
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IDENTIFICATION 

TITLE: Solution of Simultaneous Linear Algebraic Equations 
CO-OP ID: F2 UTEX LINEQN 

CATEGORY: Simultaneous Linear Equations 
PROGRAMMER: C. B. Bailey | 

DATE: August 14, 1961 


PURPOSE 
Solve one or more sets of linear algebraic equations using Gaussian 
elimination with row pivoting and back substitution. 
USAGE 
1. Calling sequence: 
The program is called by the program execute card, i.e. LINEQN.., 
in the normal sequence of Fortran control cards. 
2. Arguments: 
The following parameters and data are read inon cards. (See 9a). 
a. N - the order of the matrix in the equations Ax = b, that is, the 
number of linear equations. 
b. M - the number of vectors b for which solution vectors x are to 
be obtained, that is, the number of sets of linear equations. 
c. EP - Matrix condition parameter. (See Mathematical Method. ) 
A - the elements of the matrix of coefficients of the equations. 


e. B- the elements of the column vectors db). bo, suse bt 


3. Space required: Undetermined. 


Temporary storage required: Space is reserved for solving 50 equa- 
tions in 50 unknowns for 60 vectors b. This amounts to 8500 locations. 
5. Alarms or print-outs: If the equations are inconsistent or dependent, 
MATRIX SINGULAR is printed. 
Error returns: None. 
Error stops: None. 


Input and output tape mounting: Not applicable. 


Oo Oy as 


Input and output formats: For a more complete description of the 


I, E, and F formats, see 15. 


FZ UTEX LINEOQN - l 
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Input (data cards) 
(1) First card 


Columns 1-5 


contain the value of N ini5 format. Thus N is 
is punched a right-~-justified fixed-point integer. 


Leading zeros need not be punched. 


Columns 6-10 contain the value of M inI5 format. 


Columns 11-20 contain the value of EP in El0.4 format. 


To enter the value of 1078 one can punch 
.1E-07 right-justified in the field. 


(2) The coefficients of the matrix are read a row at atime. Each 


row is begun ona new card. Five coefficients per card are in 
16 column fields (1-16, 17-32, 33-48, 49-54, 35-80) using F16.8 


format. 


The matrix of coefficients has the following oboe | 


at 


2) 


id 


a 


a 


12 IN 
252 wee ony 
ONn2 ANN’ 


Thus, for N = 6, the input cards would contain the following 


information. | 


Card No. 1-16 


2 a 


i 

: 916 
4 Ao) 
“ 226 
6 231 
i 236 
12 2ey 
13 26 6 


Columns 

Sey 33-48 
SOU tee Scale 
aay Jy 253 
wy? "35 
262 263 


F2 UTEX LINEQN - 2 


oy O> 


49-54 55-80 
A14 415 
A204 “25 
mA PRs 
264 265 





(3) The elements of the column vectors db), bo, chee 


row wise in F 16.8 format. For M = 6, the vectors b 


» D are also read 
m 


b b 


2 ae 


are assumed to be arranged in the following manner: 


by 


Np 


ne 


_b 


boy 


522 


26 


63 


D62 


7a: 


where b. denotes the j-th component of the column vector b,. 


Thus, for N=6and M = 6, the data cards would contain the 


following information: 


Card No. 1-16 

14 by 

Ld bey 

16 bio 
17 bes 
24 Dig 
25 beg 


17-32 


boy 


22 


Do6 


Columns 
33-48 49-54 
ey Day 
bao D4 
bag O46 


55-80 


bey 


52 


Bea 


(4) Several sets of data may be processed at one run. To terminate 


the run, place one blank card after the last set of data. 


Output 


The components of the x vectors are printed row wise six per 


line in E20.11 format. 


For the above example of N=6and M = 6, 


the output values would be arranged in the following manner: 


F2 UTEX LINEQN - 3 


TOL 





EO. 
1 
‘WAR 
a3. 
14. 
LD 


Line No. Values 


J 9 Xo4 X31 0 Xqy X51 oy 
% 2 re) eos ha 52 X62 
6 X16 26 X36 X46 25 6 X66 


where Ks denotes the j-th component of the column vector x; . 


Selective jump and stop settings: Not applicable. 

Timing: Undetermined. 

Accuracy: Not applicable. 

Cautions to user: None. 

Equipment configuration: Not applicable. 

References: 

(a) Fortran System for the Control Data 1604 Computer, Control Data 
Corporation, Computer Division Publication 087A, Minneapolis, 
Minnesota (1961). 


(b) Fortran Automatic Coding System for the IBM Data Processing System, 
International Business Machines Corporation (1958). 


(c). Fortran IJ for the IBM 704 Data Processing System, International 
Business Machines Corporation (1958). 


(d) Kunz, K. S., Numerical Analysis, McGraw-Hill Book Co., Inc., 
New York, 1957. 


(e)) Fadeeva, V. N., Computational Methods of Linear Algebra, trans- 
lated by Curtis D. Benster, Dover Publications, Inc., New York, 


noo . 


MATHEMATICAL METHOD 


Gauss's method of elimination with row pivoting and back substitution is 


used, 


Suppose we wish to solve 


Ax = b and Ax=c 


for the same matrix A. This routine is designed to give the solutions for 


One or more column vectors b, c, etc. simultaneously. 


F2 UTEX LINEQN - 4 
£05 








Consider a 3 by 3 matrix A augmented by two column vectors b 


and c. ; 
Pe: “12 Zag Da Si 
See | moh oe: 223 De 25 
| eae ake 2 3 


The A matrix is triangularized by adding multiples of one row to another 
row. For column 1] find the row containing the largest element in column 1 
on or below the principal diagonal. Interchange this row with the first row. 


Subtract from row i (i= 2, 3) — times rowl. — 
ld 


_— 451 _ | 
tt aaa - cae J 1S oes 5. ; 


This eliminates the elements in column 1 below the principal diagonal. 
This process of pivoting and then eliminating elements below the diagonal 
may be repeated now on the second and third rows. The method is ex- 
tendible to any number of equations. 

Although only the A matrix is triangularized, the process of row 
interchange and addition of multiples of rows is applied to the n by ntm 
matrix consisting of the n by n coefficiént matrix augmented by m 
column vectors. 


Thus we have 
diy dio di 3 “il az 
: d52 225 “210 Maz 
: eae = 3s 
Back substitution is used. 


= 9 
ein) = i jel, ...,M 
ii i =n, .? l ‘ 
where 
S =. 0 fori=zN 
and t N 
‘Ss > di. Hej for 1°] Nel 222 


k=i+l 
The back substitution is also extendible to any order matrix. 


F2 UTEX LINEQN - 5 
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If the equations are either inconsistent or dependent (or nearly so), 
then some diagonal element of the triangularized matrix will be zero or 
nearly zero. In this case exponential overflow or even division by zero 
may occur in the back substitution phase. The alarm print-out (see C.5) 
is set off when a diagonal element of the triangularized matrix is less 
than or equal to argument EP (see C.2.c). Thus, EP is that value which 
the user is willing to call zero. It should be larger than the round-off 


error which can occur. 


FZ UTEX LINEQN - 6 
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APPENDIX B 


DETERIVINATION OF INDIVIDUAL MLMBER 
FLEXIBILITIES REQUIRED FOR BAWRIEX [F | 


Nes Devezcpment of Equations 

The symmetric matrix [F | of equation (2) is obtained 
by formulating the strair, energy in terms of internal 
forces. For example, the un-tapered bar (A, = Az) of 
Fig. Bl is attached to a web, and arted upon by an external 


axial Joad, P, and an internal axial lead, a. 


») 
g VT oar 
by 
_ a 
eer ie a4 
a "= 
1 owe 
p= 
Fie. Bl 


The strain energy of such a bar may be expressed 


(Ref. &) in terms of the end loads as 


U = |pty) dy (Bia) 
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For a load equal to 


p(y) = u + (we w) F 


the strain energy becomes 


ca Ig ue + uv + v® ) (Bib) 
345 
Upon piacing the bar of Fig. Bl in equilibrium the 
end loads may be expressed in terms of the internal forces 


as 


cS 
Hf 


y 2 30 a, + ag (B2b) 


By substituting equations (B2) into (Bl) the internal 


Strain energy can be written as 


amu = 2 ( uiea,® ~ 6 daja, + 3a,2) (B3) 
3A im 


Now referring again to equation (2), 


rs 


F F 
17 | 
F,. me {Ay oo0Ag} (2) 


re 22. 


a = [az seooas| 





it ¢an be shown by taking partial derivative of the 
energy that the IF matrix terms are merely coefficients 
of squared and ¢ress terms Buch as; 


2 BU = a Fae Fy, a, a,a,+ at F348483 + a, P3 (2a) 


Due tc the symmetry required of [F | » Fiz 1S exaceby 
one~half the coefficient of the cross-term, 21 A306 


Three factors, fyq, f45, and f55 are now introduced, 


q 


ie = f = f 
11 — 22 
3A 
= | (B4) 
r, oll 
6h. 


| 


where E45 and fo. will refer to the squared terms, and 
Pio refers to the cross-terms. 
The coefficients of a, and a3 in equations \B2) 


| are now identified by 


u = CMa, + (ka 
{ 5 } 
¥ ,* (B5) 
| v =P * 558 
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so that the member flexibilities can be found from 


; 2 
ra 
at By £1 9% ff : foo 3 
. ry ge wD eum tees 
a fh H+ toh, Re) 


= uae eo [aoe ba e 
"13 oe tooh K* fy 8, OO) 


Using the values of and from equations (B2) 


Fy ~ i 
34 
7 2 - 
Bia = = = P34 
an 
oie. 3 


II. Application to Tapered Bars 
For linearly tapered bars the coefficients £11, fy>; 
and f,, were modified by the functions of Ay/Ap, ®145 


$342, and@o>, according to the method shown in Ref. 4, 


giving: 


f? = f 
44 a4 el] 
ia ns B7 
12 12 D> (27) 
> f =x 
: 22 fo0 Boo 


fp 








Equations (B7) are then used directly in equations 
(B6). 
Consider the tapered bar shown as element number 29 


in Fig. 13, and re-drawn below in Fig. B2. The bar 


4909 | te 12 





ao ae A. = 8.36 in’ 
ee 
ve oe 
oF 
, = Cae ino Ay 
— = 463 
Ao 


Fig. B2 
geometry is obtained from Table I. Using the ratio 
Aj/Ao, the }) functions are obtained from Fig. 4 of Ref. 4, 
giving 
— Pi 44387 
3A, 
f ~ 18, > 
12 a 
t = $9» 
22 3. 
From the loading system upon the bar, as shown in Fig. 13, 


ui 


1.8001 


3.0188 


i 





the alpha (cy) and beta (/3) equations are respectively: 


ae 4.7055 ai, + 401502 a. + 7.8009 438 7 AX 


eo 








ry =< = ! 2 we Ee. ate = : 
a oth 82 he o +286 Sep O F109 


Substituting into equation (B6), Fy 3 3 becomes 
b) 


F = legals Sgt 705 500-4 4.8001 (24.7055) (= 44-83) 
1 
3313 630188) (-.4483)2 = 102,6485 


III, Applicaticn to Symmetrical Shear Panel 
For the case of the symmetrical shear panel, shown 
in Fie, B3, the strain energy can be written as in equation 


(Be): 





Fim. b3 


Symmetricai Shear Panel 


416 








| > > 
ey = | bt w, + W,) ls 78K Wy) > wad +” Goa) qi” (BS) 
l _ ! 
; ee | = ~ oy = 
Wyo re be 
ae 
and F = or U 
14 | 2 
CG G4 


A, 


pewines for elemert number 13, Fig. 13:3 


.- i aime 
or: E lr | = BE 109.7773) 
8.8 | 


IV. Application to Skewed Shear Panels 
Skewed panels are treated according to the method 


of Garvey (Ref. 10), and as used in Ref. 6. 


aT |4e5 | 
oe NB | 
& bP | 
~ 
er = 
a," { = _ 
1 24 : a 
(eo Qh 
aoe 9 mia 
Fig. BL 








The energy of the swept panel shown in Fig. Bh may be 


Written as 


ep — 2. lg ie | 
2BU = |( tt JC ay + ao) sins [1 + teat? Ti, (B9) 


Se 





where gamma iS Written as 


A Dias i 2 
j= eot™ S4ecote cot GB + cot fe 


Take for example element number 38. There is 


<2 _ 
F = et piving 
949 Cy, 


—_ 


| 2 
EF = 1 4632035.2288)(.5764)( 1 + 050505 x 2.5108)(3. 5868) 


969549 


F «= | € 78,5635) 
9,9 & 


lt should be noted that for both the bars and the 
panels, the terms cf the [F | matrix are usualiy made up 
from flexisilities of severai eiemernts. In addition, all 
main diagonal terms (F,,) are positive, All non-zero terms 
ee Ge TlexiaDreity metrix F| re given in the computer 


program data input, Avpendix C. 
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APPENDIX C 


fORTRAN PROGRAM FOR CDC 1602 
OLGITAL COMPUTER 


FORTRAN NAb:kS, SYMBOLS, AND DEFINITIONS 


A 


AT 


AR 


AT 


GJ 


BN 


BNP 


Matrix of coefficients of determinate 


Loads, basic input data, size 33x33 


ban 


nverse or Matrix A, obtained by Gauss 

3 Singjé-préecision subroutine 

Matrix of redundant loads, output, 

Size 33xl 

Matrix of all loads, output, size 49xl 
Matrix of coefficients of applied and 
redundant loads, basic input, size 49x33 
Mamnix of defleetionsdntluenee cocmiigaeats, 
Cutwnit, Sige S353 

Matrix product of (AL)(B), size 33xh9 
Negative of matrix product (H221I)(H21), 
Size 16x33 

Matrix preduct of (H22I)(H21), size 16x33 
(CJ) matrix augmented by trivial relations, 
Size 49Oxh9 

Transpose of (G) matrix 

Matrix trapie product of (GT)(UC)(G), size 
LOXLYO 


TV9 








Hil 


sur 


H2i 


UNITN 


Matrix of first 33 columms and rows of (H), 
obtained duringepartitaonine, came 3333 
Matrix of last 16 celumns (through first 

33 rows) of (H) matrix, size 33x16 

Matrix of last 16 rows (through first 33 
columns) of (H) matrix, size 16x33 
Remainder of (H) matrix, size 16x16 
Inversé of (H22) matrix 

Matrix product of (H1L2)(EN), size 33x33 
Gauss 3 subroutine error flag, where 
assigned value of 2 indicates singular matrix 
Numbering index 

Column matrix of applied loads, size 33xl 
Matrix of stress coefficients, size 49x33 
Maweix of flexi bsliies, reterred Wein 
text as (F) matrix, basic input data, 

size 49x49 

Matrix formed by placing unit diagonal 
matrix on top of matrix (EN), giving 


matrix size of 49x33 





PROGRAM BNBRGY 


* 










Clear Arrays 


Form Aa, B, UC Matrices 


Call Gauss 3 








Return 


Print 
atrix Singular 


Yes 


No 
| Invert A Matrix | 
The AI Matrix 
. | oad | 
CJ Matrix = Aq'B 


Form G Matrix From CJ By » 
Adding Trivial’ Relations 


Write Out 
A The G Matri tT 


o——+[ Form P 
Write Out 
Write Out A The H22 Matrix | 
The H Matrix | ; 
Partition H -— Form H22 


F 


u 





- “ 
* 
» = = 
° 
5 
« 


6 
a 
j 






Return 


Price 
atrix Singula 





Form :ENP 
‘| ENP = H20I x H2t 


Yes 


Call Gauss 3 


e : 
yg 


Ker « 


Invert H22 
Form H22I 


‘Write 
BN ~ENP ' \The H22T Matrix 


— 
= 
= 
<p 
* 


Form S 
41S = G x UNITN 


UNITN = I/EN’ 


Form AR 
AR = EN’x P 


a 


Write Out 
The S Matrix 

Write Out 
he AR Matrix 


O 


Form AT Write Out 
AT = S xP he AT Matrix 
° o 


Form C 
C'= 1/E(H1i1 + H12 x EN) 


Write Out 
he C Matrix. 
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APPENDIX E 


Coordinate Location of Strain Gages 


F8U-3 Wing Center Section 


Note: 1. xX = 0 = Intersection of the center section droop leading 
: edge and the center line of the airplane. 


a rs = 0 - The center line of the airplane. 


Upper Skin Inside Gages 





Gage Number | Gage Type _ x YY 
rr cP tl PN 
1-2 AX-5 63.15 5.82 
3-4 " ' 83.00 5.94 
5-6 s . 101.95 5.91 
7-8 ee T9525 5.88 
9-10 es 138. 38 5.94 : 
¥#11-12 . 153.78 5.94 
13-14-15 * AR 7-2 70.66 33.03 
16-17-18 . 94.83 33.27 
19-20-21 112.88 33.06 
22-23-24 :. 129.34 33.09 
*25-26-27 AR lor 2 153.54 33.15 
28-29-30 ie 161.02 33.33 
31-32-33 AR 7-2 127.60 75.04 
34-35-36 " ; 142.40 75.13 
37-38-39 2 155.64 75.07 
40-41-42 168.56 75.10 
43-44-45 ‘ 180.00 75.25 
46-47-48 i = 153.48 106.54 
49-50-51 ‘ F 161.20 101.26 
52253254 u ' 169.22 95.79 
55-56-57 176.49 88.92 
58-59-60 =r ” . 185.31 . 84.92 
61 A-7 144.86 * 66.10 
62 " 146.09 65.22 
63. | 7: ” 150.48 62.22 
64 ye 151.74 61.44 
65 i 146.67 68.68 
66-67-68 | AR: 722 147.90 67.87 
69-70-71 * 149.34 66.94 
72-73-74 7 150.78 65.92 


75-76-77 152.22 64.92 
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Upper Skin Inside Gages (continued) 





Gage Number — Gage Type 4 Y 


W W 
78 -A-7 153.48 64.05 
79 : 151.20 72.70 
80 " 152.61 71.74 
81 " 154.05 70.72 
82 " 155.52 | 69.76 
83-84-85 : AR 7-2 137.72 28.53 
86-87-88 Wy 140.00 32.25 
89-90-91 -." 142.55 30.48 
92-93-94 " 145.25 28.74 
95-96-97 J 147.99 45.40 
" *98-99-100 " e ' 150.30 43.84 
101-102-103 " | 152 4S" - 42.25 
104-105-106 . 155.95 58.68 
107-108-109 We, 158.11 57.21 
110-111-112 " »  ' 160.24 55.74 
113-114-115 4, 163.84 71.68 
‘116-117-118 , ot 165.77 70.51 
119-120-121 i" : 167.72 69.19 
122-123-124 " . 169.23 78.83 
“125-126-127 " 170.60 78.83 
128-129-130 " 173.00 78.83 
131-132-133 tt 147.30 110.81 
Notes 


1. a) Gages 1-2, 3-4, 5-6, 7-8, 9-10; 11-12 mounted by Chance- 
Vought and oriented with one leg. parallel to x and the 
other leg perpendicular to Yo 


b) No gage factor or other information available. 
‘2.: a) All AR 7-2 gages mounted at USNPS and oriented with one 


leg perpendicular to C.1.B., 45° gage pointing outboard - 
and toward trailing edge, other leg parallel to C.I.B. 


b) All gages from lot B-31, 120.5 + .5 ohms, 1.97 + 2%, 
b factor = -200. 
3. .a) All A-7 gages mounted at USNPS and oriented parallel to C.1.B. 
+ b) All gages from lot B-31, 120.0 + .3 ohms, 1.99 + 2%.: 
4. a) Gages 25-26-27, 28-29-30 mounted by. Chance-Vought and 
| oriented approximately the same as USNPS mounted AR 7-2 gages. 
b) No gage factor or other information available. 
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+ 


Bottom Skin and Beam Assembly 








Gage Number Gage Type X Y., 
134-135-136 AR 1 72.82 41.26 
137-138-139 a 78.08 46.75 
140-141-142 ” 92.18 : 61.29 
143-144-145 , 104.02 73.54 
146-147-148 “ia 114.54 84.50 
149-150-151 : 131.50 102.19 
152-153-154 AR 7-2 142.46 113.48 
155-156-157 AR 1 148.86 120.06 
158-159-160 AR 1 155.01 126. 39 
161-162 AX 5 70.60 O3033 
163-164-165 AR l 120.00 82.28 
166-167-168 AR 1 back-up gage to 163-165 
169-170-171 AR 7-2 146. 64 110.90 
172-173-174 AR l 91.38 40.18 
175-176-177 AR l 118.29 71.86 
178-179-180 AR 7-2 134.35 90.78 
181-182-183 AR 7-2 149.67 . 108. 86 
184-185-186 AR 1 164.74 126.64 
187-188-189 © m back-up gage to 172-174 
190-191-192 de back-up gage to 175-177 
193-194-195 4 back-up gage to 178-180 
196-197-198 AR 7-2 back-up gage to 181-183 
199-200-201 AR l back-up gage to 184-186 
202-203 AX 5 95.13 33.60 
204-205-206 AR 7-2 127.63 75.43 
207-208-209 AR 1 133.09 82.22 
210-211-212 AR 1 back-up gage to 207-209 
¥213-214-215 AR 7-2 141.50 92.64 
216-217-218 AR 7-2 152.91 106.70 
219-220-221 AR 1 110.60 42.61 
222-223-224 AR l 133.24 72.67 
225-226-227 AR 7-2 156. 76 104.11 
228-229-230 AR 1 (174.23 Tg TS 
*231-232-233 AR 1 back-up gage to 219-221 
234-235-236 AR 1 back-up gage to 222-224 
237-238-239 AR 7-2 back-up gage to 225-227 
240-241-242 AR l back-up gage to 228-230 
243-244-245 AR 7-2 113.15 33.18 
246-247-248 AR 1 146.31 82.07 
249-250-251 AR 1 back-up gage to 246-248 





Gage Number 


252-253-254 
ga 3-2 50-257 
258-259-260 


*261-262-263 


264-265-266 
267-268-269 
270-271-272 
273-274-275 
276-277-278 
279-280-281 
282-283-284 
285-286-287 
288-289-290 
291-292-293 
294-295-296 
297-298-299 
300-301-302 
303-304-305 
306-307-308 
309-310-311 
312-313-314 
315-316-317 


*318-319-320 


321-322-323 
324-325-326 
327-328-329 


*330- 331-332 


333-334-335 
336-337-338 
339-340-341 
342-343-344 
345-346-347 
348 

349-350-351 
352-353-354 
355-356-357 
358-359-360 
361-362-363 
364-365-366 
367-368-369 


Bottom Skin and Beam Assembly (continued) 


Gage Type 


AR 
AR: 
AR 
AR 

. AR 

: AR 
AR 
AR 

AR 
AR 
AR 
AR 
AR 
AR 
AR 
AR 
AR 
AR 
AR 
AR 
"AR 
AR 


* 9 5 be be 


-2 


a 
‘ 
nN 


omitted b 


4 


' xX Y 
Ww W 
160.63 101.50 
126.07 41.11 
_ +146.67 72.07 
_ 164.53 98.86 
174.23 127.57 
back-up gage to 255-257 
._ back-up: gage to 258-260 
back-up gage‘to 261-263 
back-up gage to 264-266 
» 129.19 33.78, 
154.98 74.65 
159.04 82.16 
back-up gage to 285-287 
168.71 96.19 
140.63 39.40 
-160.60 72: 52 
172.61 93.42 
192.55 126.85 
back-up gage to 294-296 
back-up: gage to 297-299 
back-up gage to 300-302. 
back-up gage to 303-305 
142.55 34.41 
147.42 34. 62 
167.90 74.74 
170.00 82.16 
back-up gage to 327-329 
177.00 90.60 
155.64 41.08 
173.60 73.93 
180.99 87.84 
201.95 —-126.43 
y numbering error ; 
back-up gage to 336-338. 
back-up gage to 339-341 
back-up gage to 342-344 
back-up gage to 345-347 
153.81 33.27 
157.59 33.15 
159.78 33.45 
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Bottom Skin and Beam Assembly (continued) 


¢ 
Ne 





Gage Number Gage Type X Me 
370-371-372 AR 1 162.58 33.06 
373-374-375 AR 1 160.90 © 30.00 
376-377-378 AR 7-2 _ 165.88 40.81 
379-380-381 as 171.08 50.93 
382-383-384 eS _ 176.64 69.04 
385-386-387 : i ' 179.49 74.20 
: 388 - 389-390 ~ AR 1 | 184.92 82.10 
: 391-392-393 * AR 1 back-up gage to 388-390 
7 394-395-396 AR 7-2 185.10 85.13 
397-398-399 AR 1 166.57 36.18 
400-401-402 | : 170.21 43.45 
. 403-404-405 tt 179.67 62.64 
406-407-408 : us ; 184.65 72.70 
/ 409-410-411 AR 7-2 | 190.93 84.71 
| 412-413-414 AR 1 208.77 120.12 
415-416-417 ee 214.08 130.93 
: 418-419-420 " — 14..17 25.34 
| 421-422-423 “ 156.06 25.09 
| , 424-425-426 " back-up gage to 418-420 
427-428-429 a. back-up gage to 421-423 
430-431-432 u 147. 39 8.56 ° 
: 433-434-435 " back-up gage to 430-432 
436-437-438 4 97.51 10.12 
439-440-441 um back-up gage to 436-438 
442-443-444 61.05 25.68 
} ® 445-446-447 " back-up gage to 442-444 
448-449-450 AR 7-2 back-up gage to 152-154 
451-452-453 " . back-up gage to 409-411 


The following listed gages are mounted externally on 
the upper and lower skins, and are back-up gages as 
indicated, or mounted over an internal beam member: 





454-455-456 AR 7-2 121.45 33.61 
457-458-459 " back-up gage to 279-281 
460-461-462 —— ? 136.59 33.61 
463-464-465 | y 151.86 33.61 
466-467-468 ue back-up gage to 367-369 
469-470-471 , 165.64 33.61 
472-473-474 " @ back-up gage to 28-30 
475-476-477 " back-up gage to 89-91 
478-479-480 - back-up gage to 19-21 
481-482-483 tt back-up gage to 13-15 
back-up gage to 161-162 


484-485-486 Jia 


wile 
Ca 
™“J 
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Bottom Skin and Beam Assembly (continued) 


Notes en 
BEAM GAGES <a 
1. a) All AR 1 gages mounted on vertical webs of beams, 


b) 


c) 


d) 


e) 


outboard of pivot rib, were mounted by Chance-Vought 
and are oriented with center leg of strain gage 
pointing inboard and on the center line of the web. 
Other legs 45° each side. (Note: Gage 178-179-180 
is an AR 7-2 mounted at USNPS, but oriented same as 


above gages.) 


All AR 1 gages mounted on vertical web of pivot rib 
were mounted by Chance-Vought and are oriented with 
center leg of strain gage pointing aft and on the center 
line of the web. Other legs 45° each side. 


All AR 1 gages mounted on vertical webs of. center 
section beams were mounted by Chance-Vought with center 
leg of strain gage pointing inboatd and on the center 


line of the web. Other legs 45° each side. 


All AR 7-2 gages mounted on vertical webs of beams, 
outboard of pivot rib, were mounted at the USNPS and are 
oriented with one leg perpendicular to CIB, and other 
legs pointing outboard and up or down depending on 
whether it is a "walk-around" gage or a "back-up" gage. 


All AR 1 gages mounted on vertical: web of Intermediate 
Rib were mounted by Chance-Vought and are oriented 
with the center leg of the strain gage pointing forward 
and on the center line of the web. Other legs 459 
each side.’ 


LOWER SKIN 


l. a) 


AX 5 gages were mounted by Chance-Vought and are 


oriented with one leg parallel to the center line of 


the skin bétween the beams and other leg rotated 
90° away. 
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Bottom Skin and Beam Assembly (continued) 


b) All AR 1 gages (361 thru 375) were mounted by Chance- 
Vought and are oriented with one leg (the first number 
on each rosette) parallel to the center line of the skin 
between the beams; the second leg (second number of 
sequence) rotated 45° counter-clockwise; and the third 
leg (third number of sequence) rotated another 45° 
counter-clockwise, making it perpendicular to first 
leg. The entire rosette points inboard. 


c) All AR 7-2 gages were mounted at USNPS and oriented with 
one leg perpendicular’ to CIB and the entire rosette 
pointing outboard and forward. 


2. No gage factor or other information available on gages mounted 
by Chance-Vought on either the beams or the lower skin. 


3. All AR 7-2. gages mounted at USNPS are from lot B-31, 
120.5 + 0.5 ohms, 1.97 + 2% gage factor, b Factor * -200. 


“The following gage elements are inoperable or give questionable 
results: 

11 
25-26-27 
98-99-100 
213 
233 
261 
318-319-320 
332 
445-446-447 


4 
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APPLBNDIX FP 
STRAIN GaGu INSTRUM. TATION 

The methed of taking strain recdinzs employed wemmon 
wheatstone bridge circuitry showm in Pie Fl. The unique 
part of the method was the use of a digital counter to 
indicate bridge unbalance and give a visual display of the 
strain directly in units of micro-inches per inch, 

For the special case considered nere where all four 
lees of the bridzse have the same known resistances R, and 
for a small change of resistance in cne leg “£4R, it can be 
showm that the standard Vheatstone bridge equation for out- 
pete yoltace Vo reduces £0 


V Ro AR 
v5 = : (F1) 


EL ES Oe BO Ne el. 


4 R (R + Ro) 


where R, is a constant resistance across terminals AC and V 
is the constant voltage source. Further, the expression 
for the gage factor (GF) that relates strain and change of 


resistance is 


(aprye = OF , (F2) 
R 


By substituting equation (F2) inte (Fl) it is immediately 
apparent that the output voltage is a linear function of 


serain, 
V (GF) Ro E 


dy. (R + Ro) 


= (Gomccantiee = (F3) 
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Utilizing this fact it was then easy to calibrate the 


electronic counter in units of strain. 


Considering again the basic bridge cireuic, assim: 


resistance R, across one of the legs effectively unbalances 


the bridge. Its effect on the basic balanced bridge 


eauations is developed below, 


(GF) € - QE 
R 


AR = R (GP)e 








on 
1,1, R+R 
R. Be 
. R Re Re 
AR=R- Ho = (F5) 


Substituting equation (F5) into (F4) and solving for Rg, 





Re = R rises = ] (F6) 


- 


From equation (F6) the shunt resistance necessary to 
produce an equivalent strain of 0.001 in. per in. when R 


is 120 ohms and (GF) equals 2.0 is then, 
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| 
Ro = 120 | —---—--—— = 4} = 59780 ohms. 
L(O0.001)(2.0) J 


If on the other hand a more standard 60000 ohm resistor 
were used for R., the equivalent strain would be 998 

micro inches per inch. This featuresweaseased tomeqlt bears 
the electronic counter. 

A Nheatstone bridge circuit was permanently constructed 
into which precision 120 ohm (2 , 25%) nesistors were 
readily connected to form the four legs. A precision 
60000 ohm (2 .05%) shunt resistor was connected with a 
Switch as shown in Fig. Fl. In addition a precision 50000 
ohm potentiometer was provided through connection CE as a 
variable resistance for any delicate balancing, 

With the switching and batancing unit and the switch F 
disconnected, the amplifier connected across terminals AC 
and the variable resistance connected acrcss CE, we have 
the special bridge circuit described previously. The 
variable resistance was then adjusted to exactly balance 
the bridge and conseauently give a zero indication on the 
electronic counter. Switch Ff was then closed and the 
amplifier gain adjusted to give a 998 reading on the counter. 
The counter had been previously set to repeatably indicate 


the integrative result of a one second sanpling of V, which 


greatly enhanced smooth readings. 
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Repeating this procedure and noting the consistancy 
of returning to zero with switch F oven and to 998 with 
the switch closed insured the establishment of a linear 
calibration of the counter to read in units of micro-inches 
per anen. 

After calibrating the counter, switches E and F were 
opened and resistances Ry and Ro were replaced by the 
switching and balancins unit leads for the active and 
compensating gages respectively as shown in Fig. Fl. The 
active and compensating gages also had 120 ohm resistances. 
Bach of the 20 pages connected through the switching and 
balancing unit was then zeroed on the switching and balancing 
unit in the conventional manner prior to subjecting them to 
load, At any time during the test run switch F could be 


closed and a change of 998 observed on the counter. This 


eave a continual check on the calibration of the counter. 
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Fig. Ft 


SCHEMATIC DIAGRAM 
STRAIN GAGE INSTRUMENTATION 


20 Channel 
Switch. & | DC Pwr: Supply 
eo bale Re . : Model 721A: 
Baldwin - Hewlett- hat Bs: 


Lima-Hamilton 











Packard. 
~ Oe 


‘ 
a in ww A 





Amplifier 
Model 111BF 


VIinvel 


Volt.-to-freq. 
Converter 
Model DY=-2210 
Dymec, Inc. 





Electronic Counter 
Model 521DL - 
Hewlett Packard 


Ry = Ro = R3 = Ry = CRCA 120 ohm (.25%) 215-RL 
Re = CRCA 60000 ohm (.05%) 215-RL 
R, = 50000 ohm Helipot TP Precision Potentiometer 


Resistance Tolerance: 
Linearity Tolerance: .5% 


145 











API iNDIX G 


WALERIMENTAL PRINCIPAL STRESSES AND 
AXIS ORIENTATION 


The magnitudes of principal stresses anda the orientation 
of the principal axes were calculated for all rosette 
locations of interest. txperimental strain readings are 
listed in Table Gl. Calculations were performed on the 
C.D.C. 1604 Digital Computer utilizing a FORTRAN program 
named ROSRED shown in Tables G2 and G3. The computer output 
is listed in Table Gil,. 

The equations solved by the program were the relations 
between rectangular rosette readings and principal stresses 


found in any standard text on the subject. 











Onan = = | oes + —__ 2€-(6, +6] ° 

max 5 fess Ces ia f 5) 

Fa = Sheen 6! (€; -€3)° 7 be, -€, +€)] ; 
2 1(1-f) (1+) 





oa & . | 
Tmax = d(1+1L) Wi (€, =€2)- + |26-€, +€3)|* 


@p = + arctan 202 = Aimer 
; C eae 


it 
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Pp is the angle from the axis of €1 to the maximum normal 
stress axis and C3, C5 and C3 are strain readings in 

the respectively numbered eages shown in Fig. Gl. Rosettes 
used in this exveriment were numbered with three cone 
secutive numbers and the lowest number identified the 
rosette. Logically then QD, is the angie from the axis of 
the lowest numbered gage in each rosette to the maximum 
normal stress axis, A positive value indicates an angle 

in the direction of €5., 

The input to program ROSRED (Table Gl) was designed 
to accommodate this notation. Only the rosette identifying 
number is listed which is also the gage number for which 
the first column of strains are listed. The second and 
third columns are then the values of strain on the next 
two consecutive numbered grapes comprizing the rosette. 

The results listed in Table G4 are identified by the 
rosette numbers which quickly orients the axis reference 
for Dp. 

It should be pointed out that the program was designed 
for specific values of E and 4, namely 10.3 x 106 psi 


and 0.32 respectively. 


14-7 












Big. G2 
FLOW CHART FOR PROGRAM ROSRED 


| Steaat 


- READ Rosette Number, 
| E149 €2i9 ©34 


a 
Rj = [B,2 + (2 Eo, - Ay)? 


G-. = 103 | As, i Rd 
0.68 1.30 7 


max 2) 





Return —no— < yes 


PRINT Rosette Number 
Omax? Omin? T max? Dp 
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Table Gi 


RIMENTAL STRAIN READINGS 


PR 


BX 


(Input To ROSRED Progran) 
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(Table Gi Continued) 
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TABLE G2 


LIST OF SYMBOLS USED IN 
ROSRED PROGRAL: 


Computer Definition 

Coded 

Name 

BPI €,, Strain in lowest numbered gage of rosette 
BP2 €,, strain in diagonal gage 

EP3 C3, strain in perpendicular gage 

STGMAX Om maximum principal stress 

SIGMIN yiny Minimum principal stress 

TAUMAX Tox maximum shearing stress 

PAIP QD , angle from the €3 axis COMO te, AX1S « 

GAGE Rosette number (lowest numbered gage in rosette) 
N Number of rosettes furnished es data input 
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Table G3 
FORTRAN PROGRAM 


"ROS RED" 
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